Modified Semi- Classical Methods for Nonlinear Quantum 

Oscillations Problems 



(N 

o 

(N 



Vincent Moncrief 
Department of Physics and Department of Mathematics, 
Yale University, P.O. Box 208120, New Haven, CT 06520, USA. 
E-mail address: vincent.moncrief@yale.edu 

Antonella Marini 

Department of Mathematics, 



■ Yeshiva University, 500 West 185th Street, 

New York, NY 10033, USA. 

■ 

and 



Department of Mathematics, 
University of L'Aquila, Via Vetoio, 
67010 L'Aquila, AQ ITALY. 
E-mail address: marini@yu.edu 

m 

»/~) \ Rachel Maitra 

■ Department of Physics, 

Albion College, 611 E. Porter Street, 
Albion, MI 49224, USA. 
s_i " E-mail address: rmaitra@albion.edu 

(Dated: January 26, 2012) 

We develop a modified semi-classical approach to the approximate solution of 
Schrodinger's equation for certain nonlinear quantum oscillations problems. In 
our approach, at lowest order, the Hamilton-Jacobi equation of the conventional 
semi-classical formalism is replaced by an inverted-potential-vanishing-energy vari- 
ant thereof. With suitable smoothness, convexity and coercivity properties imposed 
on its potential energy function, we prove, using methods drawn from the calculus 
of variations together with the (Banach space) implicit function theorem, the ex- 
istence of a global, smooth 'fundamental solution' to this equation. Higher order 
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quantum corrections thereto, for both ground and excited states, can then be com- 
puted through the integration of associated systems of linear transport equations, 
derived from Schrodinger's equation, and formal expansions for the corresponding 
energy eigenvalues obtained therefrom by imposing the natural demand for smooth- 
ness on the (successively computed) quantum corrections to the eigenfunctions. For 
the special case of linear oscillators our expansions naturally truncate, reproducing 
the well-known exact solutions for the energy eigenfunctions and eigenvalues. 

As an explicit application of our methods to computable nonlinear problems, we 
calculate a number of terms in the corresponding expansions for the one-dimensional 
anharmonic oscillators of quartic, sectic, octic, and dectic types and compare the re- 
sults obtained with those of conventional Rayleigh/Schrodinger perturbation theory. 
To the orders considered (and, conjecturally, to all orders) our eigenvalue expan- 
sions agree with those of Rayleigh/Schrodinger theory whereas our wave functions 
more accurately capture the more-rapid-than-gaussian decay known to hold for the 
exact solutions to these problems. For the quartic oscillator in particular our results 
strongly suggest that both the ground state energy eigenvalue expansion and its as- 
sociated wave function expansion are Borel summable to yield natural candidates 
for the actual exact ground state solution and its energy. 

Our techniques for proving the existence of the crucial 'fundamental solution' to 
the relevant (inverted-potential-vanishing energy) Hamilton-Jacobi equation have 
the important property of admitting interesting infinite dimensional generalizations. 
In a project paralleling the present one we shall show how this basic construction 
can be carried out for the Yang-Mills equations in Minkowski spacetime. 

PACS numbers: 02.30.Mv, 02.30.Xx, 03.65.Sq 



I. INTRODUCTION 



In this article we develop a modified semi-classical approach to the approximate solution 
of certain nonlinear quantum oscillations problems. Quantum systems of nonlinear oscilla- 
tors have of course long been studied, from the semi-classical viewpoint as well many others, 
in the mathematical physics literature. Our aim here however is not simply to provide 
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alternative proofs of already known results but instead to develop mathematical methods 
that can ultimately be applied to certain systems of quantized fields. We are pursuing that 
development in parallel to the present one and, in a companion paper, will show how some 
of the fundamental constructions of this article can, in fact, also be realized for the much 
more technically demanding case of nonabelian gauge fields Since the study of finite 
dimensional systems, however, allows one to sidestep the intricate complications of regular- 
ization and renormalization, one can push the analysis of these to a deeper level than is 
readily attainable for field theoretic problems. For this reason we have been motivated to 
carry out the present, parallel study of quantum mechanical systems from the analogous, 
modified semi-classical viewpoint. 

To some extent our work on these finite dimensional systems overlaps that already exten- 
sively developed in the microlocal analysis literature j^j but it also differs from this work in 
fundamental ways that seem crucial for our ultimate, intended application of these methods 
to quantum field theory. But even at the finite dimensional level our approach offers some 
tangible advantages over that previously developed by unifying and globalizing several of 
the fundamental steps in the analysis and by doing so in a way that, in effect, lays the 
foundation for its eventual extension to infinite dimensional systems. 

Finite systems of quantum oscillators are, of course, often considered as rough approxi- 
mations to quantized fields (through, for example, Hamiltonian lattice discretization) and 
an important class of these is characterized by Schrodinger operators of the form 

H = ~WA + mV^f) 2 + A(x\ x n ). (1.1) 

2m 2 ' 

i=l 

Here is the Laplacian on R n , = Y^=i ^ ne 'frequencies' {oJi\i € [1, . . . ,n)} are 
real and strictly positive and the potential function A : R n — > R is smooth and incorporates 
the nonlinear it ies in the sense that its Taylor expansion about the origin in R n is required 
to begin at third order, i.e., to satisfy 

dA{0,...,0) d 2 A{0,...,0) 



= Vij e [l,. ..,«]. 

We shall further require that the total potential energy function V : R" — > R defined by 



(1.2) 



n 

V(x\...,x n ) := -m + A{x\ ...,x n ) (1.3) 



i=i 
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satisfy V > on R™ and have a unique, critical point corresponding to a vanishing global 
minimum at the origin so that 

V(x) := V(x\ . . . , x n ) > 1/(0, . . . , 0) = 

(1.4) 

Vx= (x\...,x n ) eM"\(0,...,0). 

Later we shall impose a certain convexity condition on V to insure the uniqueness and 
smoothness of our basic constructions and, in the event that A(x) has indefinite sign, a 
certain coercivity condition bounding its behavior from below. Finally we shall require that 
the frequencies {wj}, characterizing the quadratic term in V(x) satisfy a convenient (but 
inessential) 'non-resonance' condition that will simplify the analysis of quantum excited 
states. 

Our approach begins by seeking a ground state wave function of the form 

£(x) = N h e- S ^' h (1.5) 

where is a normalization constant and in which S^x) is real- valued and admits a formal 
power series expansion in h which we write as 

fc2 fen 

S ft (x) ~ 5 ( o) (x) + &S (x) (x) + ^StoW + • • • + ^ (B) (x) + . . . (1.6) 

(0) 

We expand the corresponding ground state energy eigenvalue E h in the analogous way, 
writing 

(0) / (0) (0) ffl (o) fin (0) \ 

E n ~ A ^ (0) + fefa + ¥ f (2) + • • • + + . . . J (1.7) 

and substitute these ansatze into the time-independent Schrodinger equation 

* (0) (0) (0) 

F^fi = SrVb (1-8) 

requiring the latter to hold, order by order, in powers of Planck's constant. 

At leading order our formulation immediately generates the 'inverted-potential-zero- 
energy' (or 'ipze' for brevity) Hamilton- Jacobi equation, 

±-VS m ■ V% - V = 0, (1.9) 

for the function 5(o)(x). Under the convexity and coercivity hypotheses alluded to above 
(and given precisely by inequalities ( 1III.20|) and (1III.22I) below) we shall prove in Sect. 
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using methods drawn from the calculus of variations, the existence and smoothness of a 
globally defined 'fundamental' solution to Eq. (II. 9j) . The higher order 'quantum corrections' 
to 5(o) (x) (i.e., the functions 5(fe)(x) for k = 1,2,...) can then be computed through the 
integration of a set of 'transport equations' for these quantities along the integral curves 
of the gradient (semi-) flow generated by SVq) (x). The natural demand for smoothness 
of these quantum corrections will force the (heretofore, undetermined) energy coefficients 

(0) (0) (0) 

{£ (o), £(i), £(2), . . . } all to take on specific, computable values. 
Excited states will then be studied by substituting the ansatz 

$t(x) = ^(x)e-W/ R (1.10) 

into the time-independent Schrodinger equation 

* (*) (*) (*) 

HiP h = E h iP h (1. 11) 

w (*) 
and formally expanding the wave function cj) h and energy eigenvalue E h in powers of ti as 

before 

(*) (*) (*) % 2 m 

(j>h - 0(0) + + 2[0(2) + • • • > ( L12 ) 

(*) (*) / (*) o) h 2 (*) \ 

E h ~ h£ h = h ( £( ) + H£ (1) + — £ (2) + ... (1.13) 

while retaining the same 'universal' factor e _5 '^ x ^ fi determined by the ground state calcula- 
tions. 

It will prove to be convenient to reexpress the equations for an excited state in terms of 

(*) (*) 

the energy 'gap' AEn = hASn and its expansion coefficients defined via 



(*) « (0) 

A£/-( := En — £?i 



f A£ (0) + nA£ (1) + - A£ (2) + • • • + -A£ (n) + . . . J 



(1.14) 



W (*) (o) 

where A£(j) := Eu\ — £uy Under a convenient, but inessential 'non resonance' (or, more 



properly, 'non-degeneracy') condition upon the frequencies (defined precisely in (III. 441) 



below) we shall show, by a sequence of arguments given in Sections (iHBj) and (HVBl) . that 

0) 

the globally smooth solutions for 0m) can each be naturally characterized by a collection of 

w 

non- negative integers m := (m 1; . . . , m n ) and, threafter refine the notation by replacing 0( O ) 
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(m) 

with 0(0). The corresponding, lowest order energy gap coefficient will prove to be 

(m) 



A£( ) = J^m^ (1.15) 



i=l 

Cm) (m) 

so that A2£(o) = hAS^ = Y^=i will coincide with the energy gap of a collection of 

pure, harmonic oscillators, the i-th one excited to its 7?vth energy level. 

f(m) ^1 

Higher order corrections, < <j>{i)\i = 1,2,... >, to the excited state wave functions can then 
be sequentially computed through the systematic integration of an appropriate set of linear 
transport equations derived from Schrodinger's equation. The corresponding, higher order 
energy gap coefficients, < A£^\i = 1,2,... >, will each be uniquely determined by the natu- 
re 1 

ral demand for regularity of the functions < 4>{i) \ ■ A subtlety of this excited state analysis is 
that one must first develop formal power series expansions for the solutions to the relevant 
transport equations and then smoothly modify these to generate the actual, globally defined 

f(m) ^1 

smooth functions < 0(^(x)|i = 1, 2, . . . >. The formal expansions however will already suffice 

f a (ra) I 1 

to determine the energy gap coefficients < A£^\i = 1, 2, . . . > which will remain unaffected 
by the subsequent smoothing operations needed to complete the construction of the wave 
functions. 

It will become clear from our detailed analysis that all of the energy coefficients (for both 
ground and excited states) are uniquely determined from the Taylor expansion coefficients 
(about the origin in MJ 1 ) of the potential energy function V(x). Thus for non-analytic 
potential functions (many of which, of course, share the same formal Taylor expansion) the 
corresponding energy coefficients could not be expected to have better than an asymptotic 
validity as h — > 0. In this limit the wave functions become more and more sharply peaked 
about the origin and the detailed behavior of the potential energy away from the origin 
becomes accordingly less and less relevant to the determination of the energy spectrum. For 
analytic potentials on the other hand (which are of course uniquely determined by their 
convergent Taylor expansions) one expects more favorable behavior and we shall review, 
at the end of Sect. (1IVj) . the precise sense in which formal expansions, of the types given 
above in Eqs. (II.6HI.7j) for the ground states and Eqs. (II. 1 2f4TTT4|) for excited states, provide 
actual asymptotic approximations to solutions of Schrodinger's equation. The sharper error 
estimates derivable for analytic potentials will validate quantitatively the intuition sketched 
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above that such methods should indeed be more favorable in the analytic case. 

In Sect. (jV]) we shall turn from general arguments to explicit calculations by applying our 
ideas concretely to the one-dimensional examples provided by quartic, sectic, octic and dectic 
anharmonic oscillators. For each of these models we have carried out the calculation of the 
ground state wave function and its energy eigenvalue to order h 25 and compared the results 
to those of conventional Rayleigh/Schrodinger perturbation theory. In each case studied 
(and conjecturally, much more generally) we find that our eigenvalue expansions agree with 
the conventional ones and we conjecture that this agreement extends to all orders. For the 
oscillators considered it is well-known that the usual perturbative expansions for the ground 
state eigenvalues diverge but are nevertheless Borel summable to yield the exact ground 
state energies [3,4]. If the conjectured agreement of our expansions with the conventional 
ones proves valid to all orders then of course the Borel summability conclusions would apply 
to our results as well. 

For the case of the quartic oscillator in particular we present some rather striking evidence 
that the formal series (11.61) determining the ground state wave function (11.51) may itself 
be uniformly Borel summable and thus to yield a natural candidate for the exact ground 
state solution. For the quartic oscillator we also consider excited states, for an arbitrary 
excitation level, but only up through order h 3 . Here too the eigenvalue expansions agree 
precisely with those of conventional perturbation theory and again we conjecture that this 
agreement should persist to all orders. For both ground and excited states, even though 
our eigenvalue expansions agree (to the order computed at least) with the conventional 
ones, our wave functions capture, even at lowest order, the more-rapid-than-gaussian decay 
known to be valid for such anharmonic oscillators. By contrast the conventional theory 
which approximates these solutions by expansions in harmonic oscillator wave functions, 
could only hope to recapture this rapid decay in the limit that these series are somehow 
fully summed. On the other hand we suggest a procedure whereby our wave functions 
could be expanded (in the oscillator coupling constant) and truncated in such a way as to 
(conjecturally) reproduce the approximate wave functions of the conventional theory. 

To describe the ways in which our methods differ fundamentally from those of the mi- 
crolocal analysis literature will proceed more easily after we have presented our results in 
detail. We therefore postpone this comparison until the concluding section wherein we also 
discuss the (closely related) program for extending our ideas to the infinite dimensional 
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setting of quantum field theory. 



II. FINITE DIMENSIONAL OSCILLATOR SYSTEMS 



A. Ground State Preliminaries 

Substituting the ansatz (1.5) into the time-independent Schrodinger equation for the 

(0) 

ground state (1.8), expanding Sn and E% as in (1.6) and (1.7) respectively and requiring that 
the resulting formula hold order by order in h leads to the following sequence of equations 

(0) (0) (0) 

for the unknowns <J S {0) , S {1) , S( 2 ), ■ ■ ■ , £(o) , , £(2) 

^VS ( o) ■ VS ( o) 

m 2m 

--VS ( o) • vs {2) - -VS {1) ■ VS (1) + -< 
m m m 

--VS ( o) • VS m - -V£ (1) • VS (2) + ^-WA5, 
m m 2m 



V 


= o, 


(Il.la) 




(0) 




(0) 


= £(o), 


(Il.lb) 




(0) 




(i) 


= 2£ ( i), 


(II.lc) 




(0) 




(2) 


= 3£( 2 ), 


(Il.ld) 



and, for arbitrary k > 2 



fc-i 



1 1 /r' k (°) 



These consist of the 'ipze' Hamilton- Jacobi equation for previously defined (III. lap and 



a sequence of 'transport' equations (!H.lbl4lL2"l) for the higher order quantum corrections 

(0) (0) 

S(i), S( 2 ), ■ ■ ■ , £(o),£(i), ■ ■ 

As already mentioned, our strategy for solving this system will be first to show that 
Eq. (IH.lal) has a canonical, smooth, globally defined 'fundamental' solution S^y Equipped 
with this solution we shall then show that the subsequent transport equations can, sequen- 
tially, be integrated to generate the quantum corrections {S(i), 5 , ( 2 ), . . . , } to and that 
the natural demand for global regularity of these functions forces the (heretofore, unknown) 

C to) (0) (0) ^ 

energy coefficients < £( ), £(2), ■ ■ ■ f all to take on specific, computable values. 

Since S(o) will play a central role in all of our constructions (including those for the excited 
states) let us first describe, intuitively, how we propose to generate it. As already noted, 
equation (IH.lal) is the zero-energy Hamilton- Jacobi equation for the associated (inverted 
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potential) mechanics problem whose Lagrangian L ip and (conserved) energy E ip are given 
by 

I J i p (x , . . . , x , x , . . . , x ) 

n 

:=f £(^) 2 -V^V--^) (n 3) 



^(xT + n* 1 ,...,^), 



and 



2 

i=l 



Eip{x , . . . , x , x , . . . , x ) 

n 
8=1 



(II.4) 



wherein we have, for convenience, introduced the inverted potential function Vi P {x 1 , . . . , x n ) 
defined explicitly by 

V ip (x\...,x n ) = -V{x\...,x n ). (11.5) 
The Hamiltonian corresponding to Li P is given by 

ffi p (x , . . . , x , pi , . . . , p n ) 
1 n 

: =iEw+* (II . 6) 

8=1 

so that the associated (zero-energy) Hamilton- Jacob i equation (for an 'action' function 57o)) 
results from setting 

and coincides with ( III. lap . 

In view of the conditions we have imposed on V(x 1 1 . . . , x n ) (c.f. 1.3, 1.4) the graph of 
Vi P (x l , . . . , x n ) is that of a potential 'hill' (rather than 'valley' or 'well') whose summit is the 
unique, (vanishing) global maximum lying at the origin in M. n : 

V ip (x\...,x n ) < V ip (0,...,0) = 
V (x 1 , . . . , x n ) e R n \(0, . . . , 0). 
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There is a canonical family of (vanishing ip-energy) solutions to the Euler Lagrange equations 
to attempt to construct for this problem. Suppose that, for arbitrary (x 1 , . . . , a? n ) G W l 
specified at time t — 0, one could find complementary initial data (v 1 , . . . ,v >n ) G M. n such 
that the solution curve (x 1 (t), . . . ,x n (t)) determined by 

(x 1 (0),...,x"(0)) = (S )l ,..., ( x ) ") (II.9) 



and 



satisfied 



d ^,...,^) = &\...:>) (mo) 

(JjV (JjV / 



E (xHt) x n (t)- dxl(t) dxn(t) 

E w yx (t),...,x w> _ s _,...,_ s _ 



(11.11) 



and 



t lirn^ (x\t), . . . , • • • , ^) = (0, . . . , 0; 0, . . . , 0). (11.12) 

In other words suppose that, for any initial position (x 1 , . . . , x n ), one could find a (vanishing 
ip-energy) solution curve that tended asymptotically, as t \ — oo, to the potential summit 
lying at the origin in M, n . Such a solution would necessarily, in view of its vanishing ip-energy, 
also have asymptotically vanishing velocity. If such a curve existed for each (x 1 , . . . , ( x n ) G 
R n and were uniquely determined by this data then the collection of such curves would 
comprise the 'canonical' family that we are interested in. 

We shall prove below, using the direct method of the calculus of variations and additional 
convexity and coercivity hypotheses on the potential energy function V^x 1 , . . . , x n ), that 
such a solution curve does indeed exist for arbitrary (x 1 , . . . ,x >n ) G W l and that this curve is 
always uniquely determined by the given data. It will follow from the proof that the 'action' 
integral 

S(S\ . . >) : £ M L lp (,>(*), .... *•(«), - . ^) (H.13) 

converges on each such curve and that (again using the convexity hypothesis) the resulting 
function S : M n — > R is smooth (i.e., C°°) and satisfies the ipze Hamilton- Jacobi equation 
globally on W n . We shall also find that the gradient of S generates the complementary initial 
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data (v 



. . . , v n ), mentioned above, via the formula 



(o) (o) . / dS dS 





(11.14) 



and, more generally, that the canonical family of solutions curves described above corre- 
sponds to the gradient 'semi-flow' of S defined via 



where i = l,...,n, x l (0) = x . It will be natural at this point to identify S with the 
fundamental solution, SVq) , to the ipze Hamilton- Jacobi equation that we have been seeking. 

We use the term 'semi-flow' for that generated through Eq. (III. 151) since it will turn out 
that each integral curve will typically only exist on an interval of the form (— oo, e) for some 
£ > (which depends upon the curve under study). An important feature of these curves 
is that none of them (except the trivial one having (x 1 , . . . ,x n ) = (0, . . . , 0)) can achieve 
the potential energy summit at a finite time t* but instead must only approach this point 
asymptotically as t \ — oo. The well-known reason for this is that if, at any time t*, one 
had (x l (t*), . . . , x n (t*)) = (0, . . . , 0) then, by the vanishing of its ip-energy, one would also 
have (^ dx ^ \ • • • , dx - J = (0, . . . , 0) at this instant. But this is 'initial' (at time t*) data 
for the trivial solution (x l (t), . . . , x n (t)) = (0, . . . , 0) V t e R which sits at the summit for 
all t G M and thus, by uniqueness of solutions to the Euler Lagrange equations, could not 
have had the (non-trivial) initial (at time 0) conditions assumed for it. The semi-infinite 
character of the intervals of existence for these integral solution curves will play a crucial role 
in the subsequent analysis of the transport equations (III.lbHlL^l) which, along the integral 
curves of the gradient of 5(o) = <S, can now be written 



m 



dx\t) _ OS 
dt dx l 



(11.15) 




(x\t),...,x n (t)) = 
(x\t),...,x n (t)) = 




^AS {0) -8 {0) )(x 1 (t),...,x n (t)), 



,x n (t)) 



(II.16a) 



(II.16b) 



12 



and, for general k > 2, 

-(x\t),...,x n (t)) 



dS(k) ( i 



fc 

2m 



(0) 

-k£, 



j !•'•'(/) '•"(/))• 

(11.17) 



The integrals of the right hand sides of these equations, over the semi-infinite intervals 
(— oo, 0] 3 t, will prove to converge, yielding smooth, sequentially defined formulas for the 
quantities |s , (j)(x 1 , . . . ,2? n ) — 5(^(0, . . . , 0) | i = 1, 2, . . . j if and only if the energy coefli- 



(0) (0) 

cients <! £( ),£(i), . . . [> are successively defined by 



(0) 

£(o) 

(o) 
£(1) 



(0) 

£(k-i) 



_L(n) AS 



2m 



(0) 



(0,...,0), 



J_(n) AS 

2m 



(k-i) 



1 ^4 (k 



2m^j\(k-j)\ 



(0,...,0),... 
I 



(0,...,0). 



(II.18a) 
(II.18b) 

(II.18c) 



(0) 



This choice for the Su^s, together with (arbitrary) choices for the constants of integration 
{6(1) (0, . . . , 0) | i = 1,2,...} (which, however, must, for each i, be chosen independently 
of the integral curve in the canonical collection) will thus result in a well-defined sequence 
{S'(i), 5(2) . . • • • , x n ) of smooth quantum corrections to 5(0) (x 1 , . . . , x n ). Since we have 
allowed for a normalization constant, Nn, in the formula for the ground state wave function 
(1.5) one could always choose the integration constants {5(^(0, 0) | i = 1,2,...} to 
vanish without any essential loss of generality. 

Though, by construction, our canonical solution curves for the ip-problem will always exist 
on the interval (—oo, 0] they will typically only, if continued in the positive time direction, 
extend to an interval of the form (— oo, e) for some e > (where e depends upon the curve in 
question). This follows from the more-rapid-than-quadratic decay of Vi P (x 1 , . . . ,x n ) whose 
corresponding, repulsive force drives the solution curves to infinity in a finite time. For 
this reason the gradient 'flow' generated by S( ) is only a 'semi-flow' as we have already 
mentioned. 
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B. Excited State Preliminaries 



To investigate excited states we substitute the ansatz 



£ = Zfi- 3 ^ (ii.i9) 



into the (time-independent) Schrodinger equation 



- (*) h 2 , \ « « 

2m (11.20) 

W (*) v ' 



h 



W W 
and formally expand the wave function <pn and energy eigenvalue En in powers of h as before 



(*) W %(*) ft 2 M 

0fi = 0(0) + Yj0(i) + 2[0(2) + • • • (II.21a) 

(*) W /(*) ftW ft 2 W \ 

£ n := ^ = ft (£(o) + Yj^d) + 2f^(2) + ■■■)■ (II.21b) 

The function is here chosen to coincide with that previously defined for the ground state 
and thus admits a formal expansion of the same type 

S n = S m + ^S {1) + ^S {2) + .... (11.22) 

Since, however, each of the coefficients {Su) \ i = 0, 1, 2, ... } has already been computed 
at this point and does not vary with the excited state under consideration, we refrain from 
attaching a superfluous, overhead (*) thereto. 

With the ansatz (III. 191) Schrodinger's equation (III. 201) now takes the form 

|— (n) AS, - — V5-R • VS fi + V - ft£ B j 0, + h f — -M ■ V0, - ^ (n) A0* = 0. (11.23) 
But since the ground state wave function 

(0) (0) 

tPn = N h e~ s ^ h (11.24) 

satisfies 

^AS h - -J-VSr • VSr + V = ^ (11.25) 
2m 2m 

(with, by hypothesis, the same S^) we can reexpress Eq. (III. 23ft in the simplified form 

fVS h \ w h fnU W WW 



y m 



A (*) ft . . (*) ww 

J ' V0/l ~ 2m" 0ri = A ^ fi (IL26) 
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where A£ h determines the 'energy gap' 



(*) (*) 
AS h := S h 



(o) 

S n 



(*) 



ft 2 



(*) 



A£ (0) + ^ A£(!) + ^AS m + 



(11.27) 



(*) 



(with A£(j) = — between the ground and excited states. 

Substituting the foregoing expansions into Eq. (111.26j) and requiring the latter to 



hold order by order in h leads to the following sequence of equations for the unknowns 

(*) (*) (*) (*) (*) (*) 

0(o) , 0(i) > 0(2) , • • • ; A£( ) , A£(i) , AS ( 2 ) , • • 



V5, 



(o) 



m 
\ m 



(*) (*) (*) 

V0( O ) - A5 (O )0(o) = 0, 

(*) (*) (*) M (*) 

V0 ( i) - A£( O )0(i) = A£ (1) (o) + 



2m 



(n 



(*) 



} A0( O ) - 



V5 (1 



(i) 



??? 



■ V0( O ), 



(II.28a) 
(II.28b) 



(*) (*) (*) (*) (*) (*) (*) / 1 

■ V0( 2 ) - A£ (O )0(2) = A£ (2) (o) + 2A^ (1 )0 (1 ) + 



(2) 



??? 



V 0(o) - 2 



m 



(*) 

(n) A0 (1 



\ m 



(i) 

(*) 

V0(1), 



(II.28c) 



and, for arbitrary k > 1, 

VS7ni\ W W (*) , > (*) 



fc! 



AS S? ^ V ^' } 



(*) 

V0(fc_ j ) 



(11.29) 



(*) 



In view of the presence of the constant term A£( ) occurring in the 'transport operator 



£ : 



ni 



V-AS 



(0) 



11.30) 



that characterizes the left hand sides of Eqs. (IH.28fln.29l) . these transport equations require 
a somewhat different treatment from that described in the previous section for the ground 
state. Fortunately, however, there are well-developed methods for solving such systems and 
we shall review these below in Sect. (II V B p . 

The aforementioned methods begin with the development of formal power series solu- 
tions to the relevant equations and then proceed (in non-analytic cases) to 'correct' these 
formal expansions in a systematic way so as to produce actual smooth solutions. Since 
these formal expansions, however, are straightforward to develop and since they already 
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yield, by themselves, an algorithm for the computation of the energy gap coefficients 

r (*) (*> (*) ] 

< A£ (o) , A£(!) , A£(2) , ■ • • f we shall review them here, postponing the refinements needed 

r o) w w 1 

to complete the determination of the smooth wave function coefficients < 0(o) j 0(i) > 0(2) ; • • • \ 



until Sect. fUVBj) below. 

We shall prove in Sect. flHI Bj) that our fundamental solution, S , ( )(x), to the ipze 
Hamilton- Jacobi equation has the properties 

n 

5 ( o) (x) = -mJ>^) 2 + C>(|x| 3 ) (11.31) 

and 



2 

i=l 



6>^ (0) (x) = mujjxi + 0(|x| 2 ) (11.32) 
(a/, Eqs. ( 1III.68I) and (1IH.82|) respectively). Thus, writing 

1 n 

5(o) (x) = -m} j uJi(x 1 ) 2 + mS(o) (x) (11.33) 

where 



2 

i=i 



E (0 )(x) = 0(|x| J ), a,-E (0) (x) = 0(|x|') (11.34) 



(*) 



we split the operator C defined above into linear 



dx l 



and remaining, higher order terms 



i=l 

(*) (*) (*) 



dx l dx l 

i=l 



so that C = Co + £r. 

As a first step toward solving Eq. (1II.28aj) . which now reads 



let us begin with the simpler equation 



C 0(0) = 0, (11.37) 



(*) (*) 

Cof = 0. (11.38) 
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This (linear, separable) equation is easily seen to have smooth, non-trivial solutions within 

(*) 

the space of homogeneous polynomials in the n variables {x , . . . , x n } if and only if A£( ) 
has the form 

(*) (m) " 

A£( ) = Af( ) := 2 (11.39) 
i=i 

where the {rrij} are non- negative integers such that 



m 



■ =^2m t >l. (11.40) 



(m) 

A corresponding solution, /( )(x), to 



(m)(m) 

Cof(o) = (11.41) 



with 



will then be given by 



( m ) ^ — ^ . Q ( m ) 
C :=}^u l x t —-AS {0) 

i=l 

n q n 

i=l i=l 



/Jo)(x) = (x l ) mi (x 2 ) m2 . . . {x n ) m " 



:=x m 



(11.42) 



(11.43) 



and, of course, by arbitrary constant multiples thereof. 

However, the solution (1II.43j) to flll.41j) will not be unique (up to the trivial, multiplicative 



constant) unless 

( m ) * (m) . 

A£( Q ) := 2 mjUi = A£( ) = > (11.44) 
i=i i=i 

implies = rrii V % € [1, . . . , n] since, otherwise 

%)(*):= (x 1 )"* ... (a*) 1 * (11.45) 

(m) (m') 

or, more generally, arbitrary linear combinations of /(o) and /(o) would provide additional, 
independent solutions to Eq. (111.41)) . 

(m) 

While this non- degeneracy for the energy gap coefficient A£ (o) (for a fixed, chosen m = 
(mi, . . . ,m n )) is, strictly speaking, not needed for our constructions, its breakdown would 
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necessitate the extension of our power series expansions in (integral) powers of h to allow for 



half-integral powers thereo 
given in Sect. ( 1IIII) of Ref. 



_and thus somewhat complicate the analysis (c.f., the discussion 
and further references on this issue cited therein). To sidestep 
such complications (at least temporarily) we shall henceforth impose the aforementioned 
non- degeneracy requirement as a restriction upon the 'frequencies' to be considered. 
For non-degeneracy to hold, not just for some particular choice of the 'quantum numbers' 
m = {mi, . . . , m n }, but for all allowed choices, we shall need to assume that the equation 

n 

53^ = (11.46) 

i=i 

has no, nontrivial solutions for arbitrary i = {£1, ...,£„} G Z n . 

We can now proceed to construct a (formal, power series) solution to 

(m) (m) 

C f = (11.47) 



of the form 

(m) (m) 

/ = X) f- ( IL48 ) 

k=\m\ 

m) 

Knli^vi rrn r /~\ T h n i~i t-v n /~t /-\ 

homi 

variables x = {x 1 , . . . , x n }. Note that the monomials 



where each f k belongs to the space, P^ om , of homogeneous polynomials of degree k in the n 



x k :=(x 1 ) kl ...(x") k ", (11.49) 
(with k = (k±, . . . , k n ), |k| = ^ILi an d h G N U {0}V i G [1, . . . , n]), provide a basis for 



(m) 



the eigenvectors of Cq restricted to V^ om with corresponding eigenvalues given by 

(m) 



,i=i 

which, by virtue of our non-degeneracy condition, are non- vanishing whenever |k| ^ |m|. 

(m) 

Thus, in particular, £ is a bijection on V^ om provided that k 7^ |m|. 
Returning to Eq. ( III. 471) and setting 

/|L| = (x 1 )™ 1 ■ ■ ■ {x n ) mn := x m (11.51) 



one finds that, in the sense of formal power series 

(m) (m) 



1 

OO 



C /| H = J2 9k (H.52) 

fc=|m|+l 
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(m) 



where each G V\ om . Taking /| m |+i to be the (unique in V^^ 1 ) solution to 



one then finds that 



(m) (m) 

£ /|m|+l = —g\m\+l (11.53) 



(m)/(m) (m) \ ~ 

>cf/| m |+/ w+ ij = h * ( IL54 ) 



fc=|m|+2 

(m) 



for certain hk G Vf wm . Thus, taking f\ m \ + 2 to be the unique solution (in ) ^° 

(m) (m) 

£o f\m\+2 = —h\ m \ +2 , (11.55) 



one proceeds, in this way, to construct the complete formal series solution (III. 481) to (1II.47I) 
that we have been seeking. 

Clearly inhomogeneous equations can be treated in the same way provided the inhomo- 

Cm) 

geneity contains no (formal expansion) term proportional to /| m | = (x ) mi . . . (x n ) mn (the 

(m) 

kernel of Co). One thus arrives at the 
Proposition (3.4 of Ref. 

(m) 

For every formal power series g at x = there is a unique constant A( 9 ) such 
that 

(m) (m) (m) (m) 

£ / = g- Kg) f\m\ 

has a solution in the sense of formal power series. This solution is unique up to 

Cm) 

a multiple of f\ m \. 

Applying the foregoing proposition it is now clear how to solve the sequence of trans- 
port equations (1II.28I - HT.29I) in the sense of finding formal power series expansions for the 

Cm) 

{4>(k) I k = 0, 1, . . . }. Except perhaps in the case of analytic potential functions one cannot 



expect these formal series expansions to converge but we shall see later, in Sect. (HVBj) . how 



they can all be 'corrected' to yield genuine, smooth solutions to Eqs. f l!OML29|) . 



One has however, already at the formal level, succeeded to determine the energy gap 

(m) (m) (m) (m) 

coefficients {A£( ), A£(i), A£( 2 ), . . . } in the sense that, once A£( ) has been chosen, all 

(m) 

the subsequent coefficients A£(k), for k = 1,2,..., are then uniquely determined and re- 
main unmodified by the subsequent refinements needed to construct the smooth functions 

(m) 

{<pk | k = 0, 1,2, ... }. This is a straightforward consequence of the foregoing proposition 
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(m) 

and the structure of Eqs. (111.28141091 wherein each successive AS^) provides precisely the 

(m) 

constant, A(g), needed for that step in the argument. 

Even though formal expansion methods were not needed for the ground state analysis 
they could nevertheless have been used to develop formal series solutions to Eqs. (III. 1HII.2]) 

(0) (0) 

and thereby to determine the energy coefficients {£(0), £(i), • • • } defined through Eq. ( III. 181) . 
Note, furthermore, that even the nonlinear equation ( IH.lal) can be solved formally by setting 

1 n oo 

S (0 )(x) = -m^Tu^ + ^r^x) (11.56) 

i=l k=3 

where Sfc(x) G V^ om and recalling that, in view of our hypotheses for the potential energy 
function, the latter admits a formal expansion of the form 

n oo 
i=l k=3 

where Vk G V^ om . Substituting fill. 56 j) and (III. 5TI) into Eq. (IH.lal) it is straightforward 
to verify that each successive Sfc(x) is uniquely determined by an argument that closely 
parallels that sketched above for the excited states. Using the formal expansion S7o)(x) in 
place of the exact solution (x) in Eqs. ( III. lbl4lL2|) it is easy to see that these can each be 
sequentially solved (formally) by the same techniques developed for the excited states and 

(0) 

the energy coefficients {£(&) | k = 0, 1, . . . } thereby determined along with formal expansions 
{S'( / t)(x) | k — 0, 1, 2, . . . } for the actual smooth functions {S , (/ C )(x)}. 

The point of these remarks is that, for both ground and excited states, the energy and 
energy gap coefficients together with corresponding formal expansions for the wave function 
coefficients can all be computed independently of any exact construction of the solutions to 
the relevant equations. It is not difficult to see that all of these elements are determined from 
the (formal, in the non-analytic case) Taylor expansion coefficients of the potential energy 
function V(x), i.e., from the collection {C air .. A | aij G N U {0} for i G [1, . . . , n]} with 

«9l Q l\/(x) 



Cqsi,...,! 



dx 1 ai ... dx n ' 



(11.58) 



x=0 

It follows that, for non-analytic potentials (many of which share the same Taylor expansion 
coefficients but not the same energy spectrum), the computable energy and energy gap 
coefficients, for example, could not be expected to have anything more than an asymptotic 
significance in the limit as h tends to zero. For analytic potential energies though, it is known 
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that the situation is significantly more favorable and we shall review what is known about 
the precise sense in which such constructions provide approximate solutions to Schrddinger's 
equation below at the end of Sect. (1IVI) . In Sect. (jV]), which is devoted to the study of 
certain (analytic) anharmonic oscillators in one dimension, we shall discuss some evidence 
that constructions of the type under study, combined with suitable Borel resummations, 
may even yield exact results. 

We conclude this section with a brief discussion of an alternative technique for solving 
the basic equation 



(*) O) 
C 0(o) :-- 



/VS7m\ « 
' (0) » - V - AS {0) 



0(0) = (11.59) 



\ m 

through an application of the well-known Sternberg linearization theorem [5] . For a class of 
smooth vector fields including those of the type 



m 

|2> 



i=l x ' 



(with <9jE(o)(x) = 0(|x| )), Sternberg proved that, if the 'frequencies' {ui} satisfy a certain 
'non- resonance condition' (which is implied by our non-degeneracy condition ( 1II.46I) ). then 
there is a (local) diffeomorphism, 

xi — > y = /i(x), 

(11.61) 

?/ ! = /i i (i 1 ,...,i n )V«e [l,..-,n], 
defined on a neighborhood of the origin in M. n such that, in terms of the new coordinates 
{y 1 , . . . , y n } the corresponding vector field takes the purely linear form 

Y%.V = X><A (H.62) 
m z — 4 air 

8=1 y 

Using special features of the Hamilton- Jacobi (semi-) flows generated by our particular 
vector fields we shall show, in the Appendix, that Sternberg's fi can always be extended 
to a global diffeomorphism mapping IR n = {(x 1 , . . . ,x n ) \ x l G M} to a corresponding star- 
shaped domain in the space of the {y 1 , . . . , y 11 } coordinates. That the image of M. n under this 
extended fx does not normally exhaust another copy of lR n , but only a star-shaped domain 

V S 

therein, follows from the fact, already noted, that the integral curves of — typically 
only persist for semi-infinite time intervals of the type (— 00,^(7)) whereas these curves, 
expressed in the new coordinates, each have the explicit form 

yUt) = y*JO)eP* (11.63) 
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for suitable constants {(yl(0)) \ i £ [1, • • • The boundary of the star-shaped domain in 
question is thus defined by the collection of ideal endpoints with coordinates 



= ^(Oje^M. (11.64) 



In terms of the new coordinates the regular solutions to 



(*) (*) " . Q (*) (*) (*) 

£ 0(o) = 2^^iV l 0-70(0) - A£ (0 ) 0(0) = (11.65) 



are given simply by (constant multiples of) the monomials 



%)(y):=(y 1 ) mi ...(y B ) ,B - (n.66) 



with corresponding eigenvalues 

(m) 

A£( ) = 2^ (11.67) 



(m) 

as before. Note that the {0(o)(y)} are normally (i.e., for genuinely nonlinear oscillators) 
bounded by virtue of the bounded nature of the range of the y-coordinates. 

By contrast, for the case of purely linear oscillators (for which the original x-coordinates 
are already of Sternberg type but with unbounded range) the corresponding monomials, 

(m) 

0(o) (x) = (x ) rni . . . (x n ) m ", provide the highest order (in x) terms in the usual product of 
Hermite polynomials that characterizes the excited state wave functions and it is straight- 
forward to verify that the constructions of this section merely fill in the lower order terms 
of these polynomials and then terminate, thus reproducing the well-known exact results for 
linear oscillators. More precisely, one finds, by direct calculation, that all of the higher 
corrections to 

n 

Sf(x):= 5 m^^) J (11.68) 

i=l 

vanish, leaving the familiar gaussian exponential factor common to all the states, and that 
the series expansions ( 1II.21al) and ( 1II.21bl) now terminate yielding 



( AT rr ( mUJ ll\ rr ( mUJ n^n 



and 



h (pc) = N K H mi ( x l-^x l \...H mn f J-j^x" ) (11.69) 



E h = 22lmi + -)hu>i. (11.70) 

i=i ^ ' 
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III. EXISTENCE AND SMOOTHNESS OF THE FUNDAMENTAL SOLUTION 



A. Existence and Regularity of Minimizers 



A natural approach for seeking solutions to the inverted potential (ip) dynamics problem 
formulated above is to look for minimizers of the ip action functional 



lM-= / L ip (x 1 (t),... J x n (t),x 1 (t),...,x n (t)) dt 



= \2 m % { ^ {t) ^ + W * (x * (t))2 } + A (a;1(t) ' • • • ' ^ } dt 
within the affine space of curves 

£> x := {ieH l (I,R n ) |/=(-oo,0], 
7 (t)= (x 1 (t),...,x"(t)),lim 7 (t)=x 

= (x\...,x n ) Gl"}. 



(III.l) 



(III.2) 



Here H 1 ^, M. n ) is the Sobolev space of (distributional) curves in R n equipped with the norm 



IIt(0 II i? 1 



o n 



1/2 



53 + Wi a (x«(*))' 



dt 



i=l 



(III.3) 



< oo 



and x = (x 1 , . . . , x n ) is an arbitrary but fixed right endpoint lying in M. n . From the Sobolev 
embedding theorem for if s -maps [6|, |7| one has that H 1 ^, W 1 ) is continuously embedded in 



0/ T mn\ 



l e C°(I 



0/ T m>n\ 



n{-)\\L°°(I,R") = SUp 



(III.4) 



^(x*(t)) 2 <oc 



where C°(I,M. n ) is the space of continuous curves 7 : I — > lR n , and furthermore that these 
curves automatically (as a consequence of having finite H 1 norm) 'vanish at infinity' in the 
sense that 



lim |7(t)| = lim 

t— > — CXD t— > — OO 



\ i=i 



(III.5) 



0. 
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Thus the curves in D x have their (asymptotically attained) left endpoints at the origin in 
M. n which, in our setup, coincides with the unique, global maximum of the inverted potential 
function Vi P (x 1 , . . . , x n ) = — V(x 1 , . . . , x n ). 

Strictly speaking the 'curves' in if 1 (/,lR n ) are distributions but, by virtue of the Sobolev 
embedding theorem cited above, each such distribution can be represented by a continuous 
curve which (by a slight abuse of notation) we also write as 7 : / — > M. n . For this reason one 
can meaningfully speak of the values of 7(i) (as points in K n ) for any t £ I = (—00, 0] and 
also impose the right-endpoint boundary condition (included in the definition of X> x ) that 

lini7(t) =x = (x\...,x n ) GET. (III.6) 

For a more extensive discussion of the values and boundary values of curves in ff 1 (/,]R n ), 
when J is a finite domain, see Section 2.1 of Ref. Q. That we are instead working on the 
semi-infinite domain / = (—00, 0) is taken into account, as t — > —00, by the 'vanishing at 
infinity' result mentioned above. 

The first two terms in the explicit integral formula for 2^ [7] are automatically finite for 
any curve 7 G T> x since these terms comprise simply a (positive constant) multiple of the 
squared H l (I, M n ) norm of 7. The integral of the remaining term in ^[7] is also finite since, 
by assumption, the potential function A(x l , . . . , x n ) is smooth and satisfies 

V i,j G [l,...,n]. It follows that the function A(x 1 , . . . , x n )/{x, x) w , where (x, x) w := 
Y^i=\ u i( xl ) 2 y 1S bounded on bounded subsets of W 1 and thus, since the curves 7 G P x are 
each bounded, we have that 







dt A(<y(t)) 





< / dt 



oc 



A{l{t)) 



(7(f),7(*)>„ (IIL8) 



<7(*),7(*)> 

< constant \l{-)f L 2 {im < constant |t(0Ihi(j,R'»)- 

Thus the ip action integral 2^, [7] is finite for any curve 7 G P x , for any x G W 1 and, in 
view of our requirement that 

1 n 

V{x\ ...,x n ) = -J2 mu^x 1 ) 2 + A(x 1 , ...,x n )>0, (III.9) 

i=i 

^Hp[l] > for any such curve as well. Summarizing the above we have that, for any x G M. n , 
each 7 G D x can be represented by a continuous curve in C°(I, M n ) (also called 7) for which 



24 



the ip action functional is finite and bounded from below, 



l ip :V x ^R + = [0,oo), 



(111.10a) 



and for which the boundary conditions 



lim 7 (t) = (0,...,0) 



(III. 10b) 



and 



j™ 7(0 = (a: 1 ,..., 




x G R' 



n 



(III.lOc) 



are 'built in'. 

If a minimizer for X ip exists for each x e R™ we shall (tentatively) identify the sought-after 
'fundamental solution' 5 (x) by setting 



where 7x is a minimizer corresponding to the chosen boundary point x. 

Note that if a minimizing curve, corresponding to a given x, should exist but fail to 
be unique, the value of 5(0) (x) would nevertheless be well-defined. But such a failure of 
uniqueness of the minimizing curves, if it occurred, would strongly suggest a corresponding 



Hamiltonian-Jacobi equation then its gradient at x should yield the complementary, canon- 
ical 'initial data', (x, p = grad 5(o)(x)), for the corresponding, minimizing solution curve 
'lifted' up to the phase space T*R n R 2n . But the gradient of a differentiable 5( ) (x) would 
not have the multiple values needed to accommodate multiple minimizing solution curves. 
To avoid this potential difficulty we shall later impose a certain 'convexity' condition on the 
potential function V(x 1 , . . . , x n ) that will suffice to globally exclude the occurrence of multi- 
ple minimizers. Then, with some further work we shall be able to prove that 5(0) : R n — > R 
is in fact a smooth (i.e., C°°) function and that it satisfies the ipze Hamilton- Jacobi equation 
globally on R n . 

For the purpose of appealing to some standard results in the calculus of variations it 
will be convenient to work with a slight reformulation of the problem defined above. If the 



5(o) (x) = X ip [ 7x ] 



(III.ll) 



breakdown in the differentiability of the fundamental solution 5(0) (x). 

The reason for this is that if 5( )(x) were indeed a differentiable solution to the (ipze) 
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potential function A{x 1 1 . . . , x n ) were taken to vanish then it is straightforward to show that 



the curve j x G P x defined by 



7x(t) = (^V.^V"') (III. 12) 



would be the (unique) minimizer for the corresponding linear problem having boundary data 
x G W 1 . For the nonlinear problem, when ^(a; 1 , . . . ,x n ) is reinstated, one can express any 
7 G D x as 

7 = 7x + « (HI. 13) 

where m is a curve in the fixed (i.e., independent of x) linear space 

V:= \ueH 1 (I,K n ) | /= (-oo,0], limu(t) = (0, . . . , 0) 1 (111.14) 



characterized by vanishing endpoint data. For arbitrary x 6 K" we can thus seek minimizers 
hGP for the functional 

iip^c[u] :=Ii P [jx. + u] (111.15) 

rather than (equivalently) seeking minimizes 7 G P x for the original functional X ip . Writing 
u : / — > W 1 in component form, t 1 — > (u 1 ^), . . . , u n (t)) one can easily verify that 



1 n 

i=l 

f° ( 1 ™ r i 

+ y oo |2 m ^[^ (t))2+ ^ 2(M,(t))2 ] (IIL16) 

+ A(x 1 (t),...,x n (t))| dt 

where x l (t) = x l e UJit + u l (t). The first term on the right hand side of ( 1III.16I) is independent 
of u and can thus be disregarded in minimizing the action for fixed x. Terms that would be 
bilinear in x l e^ it and u l (t) are absent by virtue of the fact that 7 X (£) = (x 1 e^ 1 *, . . . ) x n e UJnt ) 
is a critical curve (in fact a minimizer) for the linearized problem. To simplify the notation 
slightly let us define 

S x [u] := X ip ,xM = ZipLTx + u] (III. 17) 

and seek minimizers for S x within this linear (Hilbert) space T>. 

Standard results for the 'direct method' in_the calculus of variations yield the following 
fundamental theorem [Theorem 1.2.5 of Ref. 



ft: 
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Let X be a reflexive Banach space and M C X a weakly (sequentially) closed 
subset. Let f : M — > R be a coercive, weakly (sequentially) lower semicontinuous 
function on M. Then m(f) = inf^M f{x) is finite and is attained at a point 
x E M; i.e., m(f) = f{x ). 

For our problem we take M = X = V, the (reflexive) Hilbert space defined above, and seek 
to verify the coercivity and (weak, sequential) lower semicontinuity of 

S X :X?^R. (III. 18) 

for any fixed x G R n . 

Coercivity of <S X on V follows if ||w(-)lifi(/,R«) —> oo always implies <S x [w] — > oo. But <S x [w] 
can be written as 

1 - 

<S X [«] = -171^2 Ui(x 1 ) 2 
i=l 



(111.19) 

+ ^ra\H-)\\m { i,^)+ I A(x 1 (t),...,x n (t)) dt 



2 

where x l (t) = x l e ulit +u l (t). Clearly a sufficient condition that ensures coercivity (but allows 
some negativity in the function A{x x , . . . , x n )) is the requirement that 

1 n 

A(x\...,x n )>--mJ2^ i ) 2 (HI.20) 

i=i 

for some constants {Aj} such that A| < u)f V i 6 {1, . . . ,n}. Simple examples where this 
inequality could be arranged to hold nontrivially include polynomials involving indefinite 
terms of less than the maximal order. We shall not attempt here to characterize the most 
general coercive choice for A but simply assume that it has been chosen to satisfy this 
condition. For Frechet differentiable functionals (such as <S X ) a useful criterion for coercivity 
is provided by Lemma 2.5.2 of Ref. Q]. 

If, as is true in our case, the functional in question is twice continuously Frechet differen- 
tiable on its domain Banach space, then a sufficient condition for its weak lower semiconti- 
nuity (derived in Lemma 2.5.1 of Ref. [6J]) is the (non-strict Vpositivity of its second Frechet 
derivative (i.e., D 2 f(x) (h,h) > in the notation of Ref. |6j) at an arbitrary point in the 



27 



domain of /. For our problem the second Frechet derivative of S x is given explicitly by 
D 2 S x [u] ■ (h,h) 

( n - 

' h\t))' + uUh\t)) 







dt {m 

' i=l 



(111.21) 



d 2 A 

dx l dxi 



(s 1 (*),...,*»(*)) h\t)h?{t) 



where u a T> and /i G T U P w X> are arbitrary elements of T> and its tangent space at u 
respectively and, as before, x l (t) = x l e UJtt + u l (t). Clearly, a sufficient condition to ensure 
the positivity of this expression is the pointwise convexity requirement 



E 



d 2 V 

dx l dxi 



w r-K'e 



> 



dx l dxi 



r")a-' 



(111.22) 



V x = (x 1 , . . . , x n ) G M n and all £ = (£\ . . . , £ n ) G M n . We shall henceforth assume that 
A(x l , . . . ,x n ) satisfies this condition as well as that needed for coercivity. It follows from 
the fundamental theorem quoted above that a minimizer for <S X always exists. 

Uniqueness of the minimizer will follow if we can show that S x (which is globally defined 
on the convex space V) is strictly convex, i.e., that for any u,v EV, u ^ v, < A < 1 =>- 



<S x (Aw + (1 - X)v) < XS^u) + (1 - X)S x {v) 



(111.23) 



(c.f., Theorem 1.1.3 of Ref. [6]). From Theorem 2.6.1 of this same reference however one 
knows that, for a continuously Frechet differentiable functional (such as 5 X : T> — » M), 
strict convexity is equivalent to strict monotonicity of the Frechet derivative which, for our 
problem, corresponds to the following (strict) inequality, 



D5 x (m) -(u-v)- DSM -(u-v)>0 



(111.24) 



holding for all u,v G T> whenever u ^ v. Written out explicitly this is equivalent to the 
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requirement that 

/o ( n 

dt lm^2{u\t)-v l (t)y 
00 I i=i 



i=l 



(0) 

7x 



(0 + «(*)) - ^ ( ( 7xW + «(*)) ) (At) - At)) 



(111.25) 



> 



whenever m-d^O. Since elements of T> vanish at both endpoints the first integral is strictly 
positive for u 7^ v and thus a sufficient condition for the strict convexity of 5 X is that 



t=i x ' 



{a 1 - P l ) > 



(111.26) 



hold for arbitrary r],a, (3 G M n . But, noting that 



<9V dV 



n f 1 H ( rSV 

''^-h + AQ + fl-AJ/J)^-^) 



n »i 

= ^(a i -/3 i ) / d\ 

i,j=i Jo 

r l n 



(111.27) 



T (77 + Aa + (1 - A)/3) (a 1 - - 



dx l dxi 



we see that this follows automatically from our previous condition on V that 



1 T 



dx l dxi 



(111.28) 



V (x 1 , . . . , rc n ) G lR n and all . . . , £ n ) G M™. Thus the foregoing sufficient condition for 
the weak lower semicontinuity of 5 X guarantees as well the uniqueness of the minimizer for 
arbitrarily chosen x G W 1 . For this reason we can write, without ambiguity, u x for the 
unique minimizer of 5 X and thus define a real-valued function 5(0) on R™ by setting 



5(0) (x) := «S x (u x ). 



(111.29) 



5(0) is, of course, our candidate for the sought-after fundamental solution to the ipze 
Hamilton- Jacobi equation but, at this point, we don't even know if 5( ) is differentiable. 
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We shall prove below that it is actually smooth and that it does indeed solve the Hamilton- 
Jacobi equation but, to lay the groundwork for this, we first need to show that the mini- 
mizing curves {ii x e £>|x e R n } are smooth solutions of the (ip) Euler Lagrange equations, 
each having vanishing (ip) energy. 

We already know that each minimizer w x can be represented by a bounded, continuous 
curve on / = (—00, 0) that, furthermore, has the asymptotic behavior 

lim m x = (0, . . . , 0) (III.30a) 

t— >— 00 

fan it* = (x\...,x n ) (III.30b) 

but, as an element of H 1 ^, IR n ), it need not, a priori, be smooth enough to satisfy the Euler 
Lagrange equations, 

dV dV ip dV 
m ^ = -^ = (IIL31) 

The machinery needed to prove that our minimizers are indeed actual smooth curves on / 
(i.e., elements of C°°(I, M n )) that further do satisfy the Euler Lagrange equations is provided, 



in technical detail, by the argument in Section 4.1 of Ref. [z|. The basic tool needed for this 
is Theorem 4.1 of the foregoing reference but, since several of the hypotheses of this theorem 
are not satisfied by our problem, we first need to show how the proof given in this reference 
can be modified to apply to the setup dealt with here. 

First of all the foregoing reference deals only with minimizers defined on bounded intervals 
of the type I = (a, b) whereas we need to handle the case for which a — > —00. But, taking 
b = 0, we can evaluate our minimizer u x at any point a G {—00, 0) and restrict the domain 
of definition of w x to the subinterval I a = (a, 0) whereon u x has the boundary values given 
by 

limu x (t) = u x (a) := y (III. 32a) 

t\a 

limw x (t)=x. (III.32b) 

It is now easy to see that this curve must be the (unique) minimizer for the corresponding 
restricted variational integral 

Zip,mb] : = / L ip (x 1 (t) 7 ...,x n (t) 7 x 1 (t) 1 ...,x n (t)) dt (111.33) 
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when the curves 7 G H 1 ((a, 0), W 1 ) are constrained to have the boundary values 

lini7(t)=y (III.34a) 

t\a 

lim 7 (t)=x. (III.34b) 

If this were not the case we could replace the segment of u x along the subinterval (a, 0) 
by the true minimizer for this segment to get a different curve w* : (—00, 0) — > M. n , still in 
H 1 (I, R n ) with the original boundary values, for which the value of 5 X is thereby decreased, 

<S X [<] <S x [u x ]. (111.35) 

But this is impossible since, by construction, u x was the actual minimizer for the original 
problem. Thus we can always 'localize' our problem to subintervals of the form treated by 
the given theorem in Ref. 7]. 

But this same theorem has hypotheses that, in our context, constrain the growth of the 
potential function V(x 1 , . . . , x n ) and its gradient VV A (x 1 , . . . , x n ) for large |x| (c.f., hypothe- 
ses (i) and (ii) of Theorem 4.1). Such constraints would normally not be satisfied by the 
potential functions we wish to consider. But since our minimizers are a priori bounded we 
are free to smoothly modify their associated potential energy functions outside of sufficiently 
large balls in MJ 1 so that the modified potentials satisfy the needed growth restrictions while 
the range of the original minimizer remains entirely within the unmodified domain. The orig- 
inal minimizing curve may now only provide a local minimizer for the modified variational 
problem (since a 'true' minimizer might need to traverse a modified portion of the domain 
of the potential function). But a local minimizer is all that is required by the remaining 
hypotheses of the cited theorem. 

In summary, with the modifications sketched above, one can apply Theorem 4.1 of Ref. [3] 
to the correspondingly modified version of our variational problem and use it to conclude 
that each of our original minimizers, w x : (—00, 0) — > M. n , is a globally smooth solution to 
the Euler Lagrange equations defined on / = (—00, 0). In fact if the potential energy should 
happen to be real analytic, instead of merely smooth, then the minimizers themselves would 
be real analytic [7] instead of merely C°° . 

Any smooth solution to the Euler Lagrange equations has conserved (ip) energy but we 
shall show momentarily that this conserved energy actually vanishes for any minimizer u x , 
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i.e., that 



E ip {x\t),... 



,...,x n (t))=0 



(111.36) 



where x i {t) = x^ 1 * + uijt) with E ip defined by Eq. ( 111.41) . 

This additional regularity of the solution curves follows from the fact that, since the 
minimizers satisfy the Euler Lagrange equations and since these curves lie in H X {I , IR n ), one 
can apply an argument similar to that leading to inequality (IIII.8j) to show that the force 
term in the Euler Lagrange equations (i.e., the right hand side of Eq. (1III.3 1|) ) has finite 
L 2 (I, M n )-norm. Thus the acceleration is square integrable over I and hence the solution 
curves actually belong to H 2 (I, MJ 1 ). From the 'vanishing- at- infinity' result cited earlier it 

(0) 

now follows that the velocity, = —^p- + ^p-, of any such minimizer vanishes as 

t — > — oo and thus that the conserved ip energy of each such curve has the value zero. 

By successively differentiating the Euler Lagrange equations with respect to t and ap- 
plying Holder's inequality and the Sobolev embedding theorem inductively to the higher 
order 'forcing terms' generated thereby, one can easily show that every minimizing solution 
curve actually belongs to H 8 (I, W 1 ) for arbitrary (integral) s > 0. Thus the acceleration 
and all the (successively computed) higher t derivatives of an arbitrary minimizer vanish as 
t — > — oo, much as for the minimizer, 



To show that 5(0) : M. n — > R (defined by flIII.291) ) is smooth (i.e., C°°) we proceed by first 
showing that the minimizing solution curves it x : / — > M n depend smoothly upon x. It will 
then follow that 




of the corresponding linearized problem. 



B. Smoothness and Asymptotics of the Fundamental Solution 



s (o)(x) ■= S x (w x ) = Zipfrx] 



(111.37) 



is defined by the composition of smooth maps and thus is smooth. 
The space 




(111.38) 




e Wnt ),s > 1 (integral) 
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is an n- dimensional, closed subspace of H s (I,M. n ) comprised of the minimizers for the lin- 
earized problem. Let 



V s :={ue H s (I,R n ) \ limw(t) = (0, . . . , 0) , 

k (111.39) 
s > 1 (integral) > 

designate the space of if ^-curves on I = (— oo, 0] having vanishing boundary data. Now, for 
any s > 3, define the Euler Lagrange map 

8 s : M s x V s -> S s " 2 , (111.40) 



where 



by 



{aGF- 2 (/,l n )}, (111.41) 



£ s ( ( i, «)*(*) := m^ffiW + u\t)) - gffUt) + u(f)) 

= m ^t^ " \ mUiU {t) + 9^ (7x(t) + U{t)) 
Using the tools discussed above it is straightforward to verify that the right hand side of 

this formula is indeed an element of X s-2 . In particular, note that, for the special case 
of polynomial potentials, the fact that the 'forcing term' f^(7x( - ) + w(-)) lies in H s (I,M. n ) 
(hence, a fortiori in H S ~ 2 (I, M n )) follows from the fact that iP-maps in one-dimension from 
an algebra under pointwise multiplication for any s > 1 [8|. Using the same tools one also 
verifies that S s is a Frechet smooth (i.e., C°°) functional of its arguments. 
From the main results of this section we know that the equation 

£•(%,«)(•) = (IH.43) 

has a unique solution, m x G T> s C T>, for any 7 X G M s Vs > 1. To show that u x depends 
smoothly on 7 X (hence smoothly on x since 7 X is linear in x) we employ the Banach space 
version of the implicit function theorem [9|. The Frechet derivative of £ s (j x ,u) with respect 
to its second argument, evaluated at u — u x , is given by 

D 2 £ s (j x ,u x ) ■ h{t) 

d 2 hHt) 2) .. . 

= m^^-mu l h{t) (nL44) 

3=1 
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we need to show that the linear operator 

D 2 £ s (t, «*) : T Ux V s -> T S S - 2 (111.45) 

defines an ismorphism of the relevant tangent spaces, T Ux T> s « P s and T £ s ~ 2 « X s ~ 2 , 
V s > 3. In other words, for any a G E s_2 we need to show that the equation 

D 2 S'(t,u x )-h{-) = <r(-) (111.46) 

admits a unique solution 

MO) e © s « t Ux p s . (HI.47) 

A straightforward method for proving this is provided again by the calculus of variations 
techniques employed above for the nonlinear problem but here specialized to the linear 
(inhomogeneous) equation given by (1HL46P - Note that Eq. (1III.46j) is the Euler Lagrange 



equation for the (inhomogeneous, 'second variation') action integral, Jx jCr [h], defined by 
J x , <T [h] = ~DlS x [u x ]-(h,h) 

/o n 
-°° i=i 

We can thus prove the needed isomorphism result by showing that j7x )CT [/t] always has a 
smooth minimizer, /i X)CT , that satisfies the Euler Lagrange equation (IIII.46I) and then showing 
that such a solution is always unique. 

Note, first of all, that the convexity requirement for V^x 1 , . . . , x n ) given by inequality 
(1III.22I) implies, when evaluated along an arbitrary minimizer, u x + 7 X , that 



1 1 £)2t/ i 

1=1 i,jr' = l i=l 



for some constants x« > > 0, i G {1, . . . ,n} and V t G J. It follows that the terms in 
J7x,ct[^] that are quadratic in h are bounded above and below by (strictly positive) constant 
multiples of the squared if 1 -norm of h. Using Holder's inequality on the term linear in h it 
is easy to show that ^F x<(T [h] is finite and bounded from below V h G V. Coercivity of J^^h] 
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follows from this same result and the weak lower semicontinuity of this functional is easily 
proven by an argument that parallels that given above for the nonlinear problem. Thus, 
applying the fundamental existence theorem cited above, it follows that a minimizer always 
exists. 

For a source of o G X s-2 , with s > 3, one can again apply the regularity arguments 
given in Ref. |7J to prove that the minimizer /i x o . so obtained actually belongs to C 2 (I,M. n ) 
and satisfies the Euler Lagrange equation (IIII.46j) in the classical sense. By appealing to 



the smoothness of V(x 1 , . . . ,x n ) one can now, upon differentiating Eq. (1III.46j) arbitrarily 
many times, easily show that this minimizer lies in C k (I,M. n ) for arbitrary k and thus is 
smooth. For an analytic potential function it would further follow, from standard results on 
ordinary differential equations, that the minimizer is analytic. Finally, from the fact that o 
and its derivatives up to order s — 2 are square integrable, it follows from the Euler Lagrange 
equation and it corresponding derivatives (again appealing to Holder's inequality and the 
Sobolev embedding results) that this minimizer belongs in fact to the Sobolev space T> s . 

There is a slight subtlety however concerning uniqueness. We really want to show that the 
Euler Lagrange equation has a unique solution lying in V s and not merely that the minimizer 
is unique. However the difference between any two, hypothetically distinct solutions hi,h 2 e 
V s would satisfy the associated homogeneous equation 

D 2 S s (t,u x )-(h 1 -h 2 )(-) = 0. (111.50) 



Using the explicit formula for D 2 S s (^,u^) ■ (hi — h 2 ), contracting Eq. flIII.50j) with (hi — 




h 2 ) and integrating over I one easily shows (using the vanishing of elements of V s at the 
endpoints to justify the integration by parts) that 

f t m)-hi(t)) 2 
+ E Sh^ {t) + Mx(t))( ^ (t) - - H(t)) | (IIL51) 

i,j=l ) 

= 

which, in view of ( 1III.49I) . clearly implies that hi = h 2 . 

Assembling the results derived above we thus conclude that 

5(o) (x) :=5 x (0=2i P [7x] (HI.52) 
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is a smooth (i.e., C°°) function, globally defined on M. n . To show that 5 , ( )(x) satisfies the 
ipze Hamilton- Jacob i equation (IH.laj) we need to calculate its gradient. Utilizing the facts 
that 7x(£) := Jx(t) + w x (t) satisfies the Euler Lagrange equation and that j K (t), together 
with its derivatives, vanish as t \ — oo, one computes that 

dS(p)(x) = — dx 



8=1 



dx 



j_J l U m E(7 i W) 2 + n7(t))|(7x(-)) 



■ (d 7 x(0) 



/0 f n 

°° 1 i=i 

n 

]Tm 7x (0)^ 



(111.53) 



+ limm^7 x Wrf7xW 

i=l 



wherein we have used the observation that d 7 X (0) = 7 X (0) = cte\ Thus, as expected, 
the gradient of at x = (x 1 , . . . ,x n ) is precisely the momentum of the solution curve 
'starting' at this point, 

^^ = m 7x (0):=^(0). (111.54) 

Since we already know that each such solution curve has vanishing ip energy it follows that 
S'(o) satisfies the ipze Hamilton- Jacobi equation, 



n 1 

Y — 

i=l 



dx { 



- V(x) = 0, 



(111.55) 



globally on M. n . These solution curves realize the gradient semi- flow determined by S*(o) in 
the sense that the solutions of 

df(t) dH n 



dt 



d Pi p=V5 (0) ( 7 (t)) 

1 dS {0) (j(t)) 



(111.56) 



m dx % m 
automatically satisfy the complementary Hamilton equation 



dpijt) 
dt 



dK 



ip 



dx 1 



(V5 (0 )(7(*)),7(*)) 



_d_ 

dx 



(111.57) 
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This follows directly from noting that 



dpijt) _y &S {0) (>y(t)) 1 ggro(7(g) 

dt ^— ' dx^dx 1 m dx 3 
i=i 

2 



a A / gg(^(7(g) \ (IH.58) 
te* ^ \ dxi J 



2m dx 

3=1 

- ^ (7(0) 

where, in the last step, we have appealed to the fact that S(o)(x) satisfies the ipze Hamilton- 
Jacobi equation (lIII.55j) . The solution curves defined by this gradient semi- flow are, by 
construction, precisely those which, followed backwards in time, tend asymptotically to the 
(unstable) equilibrium lying at the peak of the inverted potential. As such they determine 
the so-called unstable manifold of this equilibrium whereas the time-reversed solution curves 
determine the corresponding stable manifold. 

One can describe the resulting picture more geometrically in phase space by remarking 
that the equation 

p, - (111.59) 

defines a smooth, global cross section of the cotangent bundle, T*R n R™ x IR n , which 
in fact, is a Lagrangian submanifold of this canonical bundle. The last statement follows 
from the fact that the graph of the gradient of a smooth function on the base manifold, R n , 
has the maximal allowed dimension for a Lagrangian submanifold (namely n) and that the 
canonical symplectic form, 

n 

uj ■= ^dx l Adp u (111.60) 
i=i 

pulled back to this submanifold, automatically vanishes: 

i=l x ' M=l 

This Lagrangian submanifold is, of course, foliated by the solution curves (naturally lifted 
to T*R n ) defining the unstable manifold of the equilibrium at (x, p) = (0, 0) whereas the 
'time-reversed' graph defined by 

d g (Q)( x ) mifi0 v 
Vi = — ^- (111.62) 

is a complementary Lagrangian submanifold of T*R™ forliated by the solution curves of the 
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stable manifold of this same equilibrium. Note that, since 



^t( a -W) ="M (.11.63) 



2m 2 -— ' \ <9x 

i=l v 

and that V(x), by assumption, vanishes only at the origin, these two Lagrangian subman- 
ifolds intersect only at the equilibrium (x, p) = (0, 0) which can be regarded as both a 
(trivial) stable and unstable solution curve. 

We conclude this section by analyzing the leading terms in the Taylor expansion of SVq) (x) 
about the origin. This will be needed for the study of the transport equations to be developed 
in the following section. From the definition of S(o) (x) we clearly have that 

KmS( O )(x) = (111.64) 

and, from the argument of the preceding paragraph, we also know that 

lim ™®H = 0. (111.65) 

Thus the Taylor expansion of 5(o) (x) about the origin begins with the second order term. 
To calculate this term precisely recall that (x) can be expressed as 

S(p) (x) = <S x (w x ) = X iPiX [w x ] (111.66) 

which, in view of Eq. ( 1111.161) results in 

1 - 

5 (o)( x ) = 2 m ^2 UJ ^ 2 
i=i 

J ^ { 2 m E {«W) 2 + (IIL67) 

+ A ^7 x (t) + Mx (t)) | dt. 

We shall find that the Taylor expansion of the integral in this formula begins at third order 
and thus that 

1 n 

5 (0) (x) = -mJ2^) 2 + 0(M 3 ). (III.68) 

i=i 

To see this we compute the gradient of expression flIII.671) directly, exploiting the facts 
that w x (t) satisfies the Euler-Lagrange equation (1III.31I) . vanishes at the upper endpoint of 
the domain of integration and has time derivative, m x (£), vanishing at the lower endpoint, 

<(0) = 0, lim w x (t) = 0. (111.69) 

t\— oo 



+ 



38 



The computation gives 



f)A 



oo 



(111.70) 



and it will follow that the final term vanishes provided that cfyw x (£) remains bounded as 
t — > — oo. But our implicit function theorem argument showed that the Euler-Lagrange 
equation (1III.31I) (implicitly) determined, for each s > 1, a smooth functional 



X s : M s ->■ X? s 



(III. 71) 



such that, reverting to our previous notation, 



u* = x (7xJ 



(III. 72) 



It follows that 



where {djj x (t)) = (0 



0). Since 9j7 x clearly lies in 



> • • • i 



(III. 73) 



T (0) M S = H s (I,M. n ) 



7x 



and since the linear operator Dx 5 ^) yields an isomorphism from this space to 



(III. 74) 



T (o) V S = T U V S ^V S , 

x s Ctx) 

it follows that clit x e and hence is not only bounded but, in fact, satisfies 



(III. 75) 



lim d jU i(t) = 0. 

t\— oo 



(III. 76) 



More explicitly djU x is the (unique, smooth) solution to the linearlized Euler-Lagrange 
equation, 



d 2 A xo) 



E^^(w)+«x(0)^) 



k=l 



dx l dx k 
d 2 A 



(III. 77) 



fc=i 
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which is guaranteed to exist by the analysis given previously. 

Thus we conclude that the expression for the gradient of S'(o) simplifies to 



dt 



dA 



(111.78) 



Differentiating again one obtains 



+ 



dt 



k=l 



dx £ 



(III. 79) 



Recalling that the Taylor expansion of A{x 1 1 . . . , x n ) about the origin begins at third order 
(c.f., Eq. (IIII.7j) ) and noting that 



x— >0 



we thus get, upon taking the indicated limit, that 



(111.80) 



lim dedjS^o)^) = mujd J e 



(111.81) 



from which Eq. 68[) then follows. In the same way one sees that the Taylor expansion 
of the integral in Eq. (1III.70j) begins at second order and thus that 



<9jS(o)(x) = mujX 3 + 0(|x| 



(111.82) 



The compatibility of these expansions with the ipze Hamilton- Jacobi equation is easily 
verified. 



IV. INTEGRATION OF THE TRANSPORT EQUATIONS 

A. Ground State Analysis 

As discussed at the end of Sect. (JTTJ) one can, for convenience, choose additive constants 
of integration in such a way that the quantum corrections, {^(x) | i — 1, 2, . . . }, to S(o)(x) 
all vanish at the origin. Making this choice one finds, from Eqs. (III. 161) and (1II.17j) that the 
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{S^)(x 1 , . . . , x n )} must be given by 



/ 1 (°) 



...,x n )= I dt ( ^ (n) A% - 5 (0) ) (7x(*) ) (IV.la) 



.Wx 1 ,...,^) 



'(2) 

and, for arbitrary /c > 2, by 
^(x 1 ,. . .,x 



f° dt (- {n) AS {1) - -VS {1) ■ V% - 2£ {1) ) ( 7x (t)) (IV.lb) 



/° f k 1 ^— i A-' (°) \ 

(IV.2) 

where, recalling Eqs. fllll. 124 - ffTTTT4|) 

7x(*) = 7x(*) + «x(0 

(IV.3) 

= (x 1 e^,...,x"e^) + KW,---,<W) 
provided that the integrals all converge. Since these integrals extend over semi-infinite ranges 
along curves that asymptotically approach the origin as t \ — oo convergence is possible only 

(0) (0) (0) 

if the energy coefficients {£(0) , £(i) > £(2) , • • • } are sequentially chosen so that the integrands 
above all vanish as x — > (i.e., according to Eqs. flll.18aflll.18cp ). In particular, in view 
of Eq. (1III.81j) , we are forced to choose 

i=i 

that the integrand f^^AS^ — £(0)^ has Taylor expansion (about the origin) beginning 

(0) 

with the first order term. If Sn) (x) exists and is smooth one can compute £^ via Eq. (III.18bj) 
and attempt to evaluate 5"( 2 )(x) through (llV.lbj) . If 5*(2)(x) exists and is smooth one can 

(0) 

calculate £(2) and proceed accordingly. 

Clearly this inductive procedure will work provided that integrals of the form 
f- dt Q ( 7x (^)) exist and yield smooth functions of x = (a; 1 , . . . , x n ) whenever £?(x) is 
smooth and satisfies Q{0) = 0. To establish that this is true we shall first need an estimate 
for the asymptotic behaviors of the curves 7 x (t) as t \ — 00. 

Recalling that the curves { 7x (£) | t G /, x 6 IR n } are precisely the integral curves of the 
gradient (semi-) flow of 5(o) defined by 

d^M _ 1 as (0) ( 7 x(t)) 

dt m dx { (IV. 5) 

7*(0) = x i , i = l,...,n 



so 
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one computes, using the ipze Hamilton- Jacobi equation satisfied by 5( )(x), that, along this 
flow, 5(o) obeys the evolution equation 

dS<f» 1 A dS {0) ( 7x (t)) dS m ( 7x (t)) 

dt 17x1 )] dx i dx i 

i=l 

= 2V( lx (t)) (IV.6) 

n 

= m^^ 2 (7x(t)) 2 + 2A( 7x (t)). 

i=l 

On the other hand we know, from Eq. ( 1III.68P and the definition of A(x), that 

1 n 

5 ( o) (x) = -m^w^x^ + Odxl 3 ) (IV.7a) 

and 



2 

i=l 



A(x) = 0(|x| 3 ) (IV. 7b) 

It follows that, on a sufficiently small ball centered at the origin and having (Euclidean) 
radius e, we have 

(7x(t)) > 2(a; min - Ce)5 (0) ( 7x (t)) (IV.8) 

for some constant C > and with w min = min{o;i, . . . , u n } > 0. Thus, choosing e sufficiently 
small to ensure that u min — Ce > 0, we get that 

5(o) (7x(t)) < K x e 2 ^- C ^ (IV.9a) 

with 

K x = 5(o) ( 7 x(£)) e - 2 ^^- C£ )*S (IV.9b) 

for all i in the range — oo < t < £* < where t* is a time prior to which 7x (i) lies entirely 
within the chosen ball of radius e. In view of Eq. (lIV.7a|) and the fact that e can be chosen 
arbitrarily small it follows from (IIV.9j) that, as t \ — oo, the Euclidean length, | 7x (£)|, of 
7x(0 decays at least exponentially rapidly in the sense that 

|7x(*)| < |7x(t x )|e (Wmi ^ Ce)( *-* ;) (IV.10) 

for all t < tt < 0. 
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With this exponentially decaying bound on |7 x (t)| it is straightforward to verify the 
convergence of integrals of the form Q(^y x (t)) dt for arbitrary smooth Q{x) that vanishes 
at x = 0. To complete the sequential construction of the {S(i)(x) | i = 1,2,...} it suffices 
to show that such integrals are automatically smooth in x. 

For this purpose it is convenient to reexpress these integrals in terms of the Sternberg 
coordinates for the vector field — V-S'(o) introduced in Sect. (Ill Bp and extended to a global 



diffeomorphism in the Appendix. From Eqs. (111.61j) and (111.63j) it is clear that 



G fa(t)) = G ° ^(y)] {t) (iv.li) 



where 



y;(t) = //(*\...,*")e^ 

(IV.12) 

= y i (x 1 ,...x n )e^* 

in which the {// | % — 1, . . . ,n} are globally smooth functions on M. n and where Q o /i _1 (y) 
is smooth throughout the (star-shaped) domain of definition of the {y % \ i — 1, . . . , n}. 
It follows that 



y=y T (*) 9x 



e^dt (IV. 13) 



Since each ojj > and since the factor multiplying each e u i is smooth and bounded along 
7 X this integral clearly converges for all x e IR n . Higher derivatives are readily computed in 
the same way and it thus follows that f° Q ("f x (t)) dt is globally smooth. 



B. Excited State Analysis 



In Section (III Bj) we showed how to construct formal Taylor expansions for the unknown, 

(m) 

excited state wave functions {(/>ik) I k — 0, 1, 2, . . .}. By exploiting well-known techniques for 
constructing globally smooth functions on M. n that have arbitrarily specified Taylor expan- 



sions about the origin 



10] one can generate a set of functions {v\k) \ k — 0, 1, 2, . . .}, smooth 



(m) 



on W 1 , that have the Taylor expansions needed for the unknowns {</>(&) | k = 0, 1, 2, . . .}. It 
follows that one can now seek solutions to Eqs. ( III.28l4n.29j) of the form 

(m) 



(m) (m) 

(fc) = V(k) + X(k) 



(IV. 14) 
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where the (x\k) I k = 0, 1,2, . . .} are smooth but vanish to infinite order (i.e., have trivial 
Taylor expansions) at the origin. Such functions are often referred to as 'flat'. 

Substituting the decomposition flIV.141) into Eqs. f llX28lhT29l) (with (*) — > (m)) and 
exploiting the fact that the by construction, formally satisfy these equations to all 

orders at the origin, one finds that the \X(k) I k = 0, 1, 2, . . .} must satisfy equations of the 
form 



fVS(Q)\ fm) (™) (m) (m) 



(IV.15) 



where the 'source' terms $j\k) \ k = 0,1,2,...} are smooth and vanish to infinite order at 



the origin. Integrating Eq. (1IV.15I) along the integral curves 7 X of the vector field ^ one 
easily finds that the only functions that could have the desired properties would be given by 



Cm) , . 

X(fc)(x) 



(m) 



(IV. 16) 



provided that the integral expressions converge to smooth functions that vanish to infinite 
order. From the exponential decay of 7 X (£) as t \ — oo (c./., inequality fllV.lOj) ) and the fact 
that each ( <T(fc)(x 1 , . . . ,x n ) vanishes to infinite order it follows that ( <7(fc) (7x(£)) = 0(e~ c '*') for 
every C > and thus that the integrals above always converge. To verify their smoothness 
and flatness properties it is convenient again to reexpress the curves 7 X in terms of Sternberg 

(m) 



coordinates. Writing xttyfe) m the equivalent form 



X(fc)( X ) 



y=y7(*). 



dt 



(IV. 17) 



where 



yUt) = ^{x 1 



y [x 



x n )e Uit 



n\ uiit 

x e 



(IV. 18) 



and exploiting the fact that Sternberg coordinates, by construction, satisfy 



X ' 



x l + 0(|x| 



(IV.19) 



it is straightforward to verify that the above expression for X(&)(x) vanishes in the limit 
x — > 0. Partial derivatives of the X(fc)( x ) can be successively computed as in the previous 
subsection and verified to vanish as above. It follows that the x|fc)(x) are smooth and vanish 
to infinite order as desired. 
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The above argument is essentially equivalent to that given in Sect. (3) of Ref. j^j]. Our 
equation for thexjfc) is slightly simpler than the one dealt with therein however since we have 
formulated the excited state problem differently by arranging to solve for the energy gaps 
rather than the actual energies. Our parametrization of the excited state wave functions 
also differs from that of Ref. Q in that we have included all the (ground state) quantum 
corrections to SVq) in the exponential factor appearing in our excited state ansatz. 



C. Asymptotic Character of the Formal Expansions 

The principal difference between our approach and that of the microlocal analysis liter- 
ature lies in the different methods used to determine the 'fundamental solution' 5 , ( )(x) to 
the ipze Hamilton- Jacobi equation. But once this function has been obtained the further 
computation (in our approach) of its quantum corrections {S(j)(x) | i — 1, 2, . . .} and, more 

(m) 

generally, the excited state wave functions{0(fc)(x) | k = 0, 1, 2, . . .} proceeds more or less 
in parallel to that discussed in [2J, allowance being made for the different parameterizations 
adopted. 

Because of the different parameterizations employed there is no direct, one-to-one corre- 



spondence between our quantum corrections and those of [2J but each set clearly determines 
the other with the same information being simply 'packaged' in different ways. Thus the 
asymptotic nature of these earlier expansions, established in {2] and other references cited 
therein, applies equally well to our expansions. For that reason we shall simply quote these 
earlier results as giving a precise meaning to the sense in which Schrodinger's equations is 
being approximately solved. 

In Ref. (2! both ground and excited states are sought in the form 

V>(x) = ( a(x,n)e- 5 (°) (x)//l (IV.20) 

and formal expansions (in powers of H) are developed for the coefficient wave functions 
( a(x, h). But then one can 'realize' these formal expansions by actual smooth functions to 
which the expansions represent asymptotic approximations, 

( a(x, K) ~ ( aj )(x) + 7w(i)(x) + ^a( 2 )(x) + . . . , (IV.21) 
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in the sense that 



<9x ||a- ^a 0) (x) fi?J 



<C KaN h N+1 , xeX 



(IV.22) 



for all K CC R n , a G N n , JVeN. Realizing, in the analogous way, the formal expansions 
for the energy levels via 



(m) (m) (m) ^2 {m) 

E(h) ~ £ (0 ) + + 2, E { 2) + ■■■ 



(IV.23) 



one proves (c.f., Theorem 3.6 of |2j) that 



(IV.24) 



where |<9£r(x, < C K , N>a h N ', x G K for every X CC M n , a e N n , N e N. For analytic 
potentials the bound that the error vanishes faster than any integral power of h can be 
replaced by the statement that it vanishes faster than e~ € °^ h for some eo > 11]. 

There are elements of arbitrariness in the procedure of passing from formal series to their 
'realizations' alluded to above so one could hardly hope to obtain actual exact solutions 
to Schrodinger's equation in this way. In the section below however we shall investigate 
some computable examples of our approach for which, through Borel resummation, one may 
perhaps be able to adduce exact results. 



V. COMPUTABLE ONE-DIMENSIONAL EXAMPLES 



A. Ground State Calculations 



As explicitly computable examples of the above framework we here study the family of 
one-dimensional anharmonic oscillators characterized by potential energies of the form 

1 



V K (x) = ^mu 2 Q x 2 + gx 2K , 



k = 2,3,4... 



(V.l) 



concentrating, for simplicity, on the quartic case (with k — 2). The ipze Hamilton- Jacobi 
equation for the quartic oscillator is readily integrated to yield 



2, ,3 




m oj, 
Qg 



1 + 



2gx' 



2\ 3/2 



mtu. 



1 



(V.2) 
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and thus to give a zeroth-order approximation to the ground state wave function of the form 



(0) 



# ^ N e - {s ^ /h) + - (V.3) 

where N is a normalization constant. This approximation already begins to capture the 
more-rapid-than-gaussian decay at large |x| expected for the actual ground state solution 
and reduces to the exact (gaussian) harmonic oscillator result in the limit g ~\ 0. 
The transport equation for Sm takes the form 

1 dS {0) dS w _ (<g _ 1 dj^ (v 4) 



m dx dx V 2m dx 2 

dS 

which, in view of the linear behavior of in x for small x, 



dSro) / 2gx 2 

— — = mu) x\ H (V.5) 

(0) 

generates a logrithmically singular solution for unless one fixes E(q) by demanding that 



£ ^ = H2^^ = r°- (v - 6) 



(0) 




Substituting this choice for £( Q ) into Eq. flV.4j) and integrating yields the (everywhere smooth) 
result 

1 f / + 

(V.7) 

No- 
where we have adjusted the arbitrary additive constant of integration to arrange that 
lim^o 5(i) =0. 

(0) s 

Thus the first order correction to ip, using ^ = + + . . . becomes 

3/2 



■2, ,3 
^0 



1 | 2gx* 
"^0 



to) N . e 69 

V> = ~ ^ : ~, 1/2 ( V - 8 ) 

1/4 / 1+ 1- 



1 + 



2gx 2 



and, at this order, the approximation for the ground state energy is simply 

(0) 1 

E n = -tkj + ■■■ (V.9) 
which of course is just the exact result for a harmonic oscillator. 
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Continuing in this way one finds the transport equation for S( 2 ) to be 



1 dS(o) dS(2) 
m dx dx 



(0) 1 

2£(i) + 



m 



dS, 



(i> 



dx 



1 d 2 S t 



(i) 



(V.10) 



for which, to avoid a logarithmic singularity in the solution, one is forced to take 



(0) 

£(i) 



1 d 2 S a 



(i) 




lim . , 

i->o 2m aar 

4m 2 Ug 

The corresponding correction for the ground state energy thus takes the form 



(V.ll) 



E h 



■ 1 3 gft 

«W0^ + - o 5 + • • • 



2 ' 4 m 2 Wn 



(V.12) 



and, upon integration, Eq. ( IV.lOj) then yields 



S 



(2) 



1 - a/1 + 



2g£ 



+ 20(7mx 2 + i^| 



2g^ 

mujf, 



3/2 



(V.13) 



6x 2 (mu;o) 3 (^1 + 

In spite of the factor 5(2) is readily verified to be everywhere smooth. 

This pattern continues indefinitely in the sense that only differentiations and eval- 
uations of smooth functions at x = are needed to calculate the energy coefficients 

(0) (0) (0) ^ 

£(o), £(i), £(2), • • • f while elementary integrations suffice to generate the sequence of quan- 
tum corrections {5(i), 5(2), 5(3), • • • } to 5(0). In particular the transport equation at order 
k yields an explicit smooth expression for in terms of previously computed quantities 
and the corresponding energy coefficient is then given by 

1 d 2 S( k ) 



(») 
£ 



(fc) 



lim 



(V.14) 



2m dx 2 

We have carried out these calculations explicitly to compute the first 26 terms in each of 
the expansions 

(0) / (0) (0) h 2 (0) h 3 (0) \ 

(V.15) 



((» 



(o) (0) fp(o) ^3 (0 ) 

h [ £ (0) + h£ {1) + —£ ( 2) + —£( 3 ) + 



and 



S R 



h 2 h 3 
5(0) + HSm + — S { 



St*) + . . . 



(V.16) 



'(0) I" T -y O {2) T — 0(3) 

and find that the energy expansion (IV. 151) agrees precisely with the conventional result 
obtained from (Rayleigh-Schrodinger) perturbation theory carried out to the corresponding 
order (c.f., [12|). 



48 



As discussed more fully below we find analogous, precise agreement with the conventional 
ground state energy expansions (carried out to the same order) for the sectic, octic and dectic 
oscillators characterized by 

V K {x) = l -mulx 2 + gx 2 * (V.17) 

with k = 3, 4, 5. In view of these results it seems plausible that agreement of the energy 
expansions extends to all orders in h for each of the oscillators studied and perhaps also 
to a much larger family of anharmonic oscillators (e.g., those having potential energies V K 
with k — 6, 7, ... or, more generally, those allowing intermediate terms of the form rx a with 
2 < a < 2k). One might hope to verify these conjectures by deriving suitable recurrence 
relations for the energy coefficients in the various cases under study and comparing them with 
correspon ding recurrence relations derivable from the conventional (Rayleigh-Schrodinger) 
approach [13| but we shall not pursue that issue here. 

For the quartic oscillator in particular the energy expansion can be rewritten as 

2 



E 



quartic 




(V.18) 



where J^uartic is a formal power series in the dimensionless quantity . It is well-known 
that this series, derived from conventional perturbation theory as a power series in the 
coupling constant g (and thus the conjecturally equivalent series derived as above as a 
power series in h), is divergent but asymptotic and amenable to Borel resummation. In 
particular the Borel sum has been proven to equal the exact ground state energy for this 
oscillator ^| and corresponding results are known to hold for the higher order oscillators as 
well [3J. The exact expression for this energy (in say the quartic case) is not analytic in 
fi := about jj, = and thus a formal series expansion for this quantity in powers of fi 
can never converge. Fortunately the asymptotic (and in fact Borel summable) character of 
this series allows one to extract accurate results from it nevertheless. 

For the sectic oscillator (k — 3 in the potential expression) one derives, as above, the 
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formal series 



(°) f 1 15 ( gfi 

-^sectic — n^O { 77 + 



2 8 U 3 w 4 
195 ( (jlr V s 1239675 / gh 



64 \m 3 ojQ J 256 ym 3 ^ 

whereas, for the octic and dectic oscillators one gets 

(o) r i 105 #/i 3 

-^octic — "^0 <( - + 



(V.19) 



2 16 \m 4 w 5 
67515 ( gh 3 V + 401548875 ( gh 3 



32 Km^ulJ 128 \m 4 a;g 



(V.20) 



and 



(o) r i 945 f #n 

-^dectic — lUJo \ - + 



2 32 V^ 5 ^o 
140057505 f gh A \ 2 78210463124745 / gh 



+ 



1024 Vm 5 a;gy 16384 \m 5 u$ 



(V.21) 



respectively. As mentioned above we have computed each of these expansions, by the method 
developed herein, up to the term of order h 25 and find that they each coincide with the 
corresponding (truncated) series expansions derived via Rayleigh-Schrodinger perterbation 
theory. Notice however that since fi sectic := J^-, fi octic : = and fi dcctic : = 4^ are the 
relevant, dimensionless expansion parameters for each of the indicated cases, and that only 
integral powers of g occur in each of the formal expansions, there are actually fewer (non- 
vanishing) terms in each successive expansion. More precisely those energy coefficients that 
would correspond to non-integral powers of g were explicitly found to vanish even though 
the corresponding quantum corrections to S<q\ were non-vanishing. 

Remarkably all of the integrals involved in computing the quantum corrections 
{S'(i), 5(2), 5(3), . . . } to 5(0) (up to the highest order computed, namely 5(25)) were express- 
ible explicitly in terms of elementary functions for the quartic and sectic oscillators whereas 
for the octic and dectic cases some (but not all) of the quantum corrections required, in 
addition, hypergeometric functions for their evaluation. It seems plausible to conjecture 
that these patterns persist to all orders in h and thus, for the quartic and sectic cases in 
particular, lead to formal expansions for 5^ in terms of elementary functions. 

The quartic case is especially interesting in that each of the computed functions 
{ ~faP '■■■■> dS <il 5) } ^ s seen ) by inspection, to be odd in x and to have uniformly definite 
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sign on each of the intervals (— oo, 0) and (0, oo). The corresponding integrals for 5( ) and 
5(i) have already been given in Eqs. (IV. 2 j) and (IV. 7j) respectively and these functions behave 
exceptionally whereas, for the subsequent quantum corrections {S(e)(x),£ = 2, 3, 4, ... , 25} 
one finds that their constants of integrations can always be chosen so that each of the func- 
tions Sm : R — > M has definite sign (alternating with £) throughout its domain and each 
decays monotonically from its global maximum (or minimum, depending upon £) at x = 
to the value as \x\ — > oo. It seems plausible to conjecture that this pattern continues 
to all orders in h for the quartic case and, if so, that one might establish such a result 
mathematically by a suitable inductive argument. 

One further finds that the sequence of ratios < — S( ( e ] — , i — 2, 3, . . . , 25 > decreases mono- 

l %) J 
tonically with increasing i. If this pattern persists and if the aforementioned conjectures for 

the quartic case can be verified it would follow not only that the formal series X]^2 (^ S{e) ^ h j 

is Borel summable but also that Y2^2 { ^ (l)<y ^ h ^ is uniformly Borel summable for all x G M. 

Under these circumstances the function 

(°) S (0)M „ , , s(x) 

ijj(x) := Nen 5 d)W— (V.22) 

where S(x) is the Borel sum of the formal series is a normalization 

constant, would be a natural candidate for the exact ground state wave function for the 
quartic oscillator. 

For the sectic and higher order oscillators studied explicitly the corresponding patterns 
are more difficult to identify but one should keep in mind that the observed and conjectured 
simple features of the quartic oscillator were sufficient but not strictly necessary for its Borel 
summability. 

B. Excited States for Quartic Oscillators 

We now turn to the construction of formal expansions for excited state wave functions, 
concentrating, for simplicity, on the quartic case. Extending our notation 

(0) (0) /(0) (0) ^2 ( o) ^3 ( o) \ 

E n := hS n = h (£ (0) + h£ (l) + — 8 {2) + — £ (3) + . . . J (V.23) 
for the ground state energy to that for a generic excited state we write 

(*) (*) /(*) {*) ft 2 M ft 3 M \ 

E h := h£ h = h \S (0) + h£ (1) + — S (2) + —£ {3) + ...). (V.24) 
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We shall soon find however that these states can, as in the case of a pure harmonic oscillator, 
be naturally labeled by a positive integer n and thenceforth sharpen and unify the above 
notation to 

(n) (n) /(«) (n) ft 2 (n) fft (n) \ 

E n := HS n = h i£ m + h£ (1) + — £ (2) + — £ (3) + . . . J (V.25) 
where n = 0, 1, 2, 3, . . .. 

(*) 

We look for excited states ^(x) by setting 

Mx) = Mx)e- Sh(x)/ \ (V.26) 

where 

S n (x) = (s (0) (x) + hS {1) (x) + + ^ (3) (x) + . . (V.27) 

is the formal expansion defined for the ground state wave function in the preceeding sec- 

(*) 

tion, and deriving the associated Schrodinger equation for the factor (j)n(x). Since Sn(x) 
will remain fixed throughout we refrain from attaching a superfluous overhead 'naught' or 

'star' to it. We shall, of course, only be interested in generating non-constant solutions to 

(*) 

the Schrodinger equation for (pn(x) since the trivial, constant solutions merely reproduce 
multiples of the ground state. 

Letting H designate the Schrodinger operator for the quartic oscillator, 

"-^ +1 2^ 2+ ^ < V ' 28 > 

we thus seek solutions to 

H U h e- Sh/n ) = f&l ($ne- Sh/ ^ (V.29) 

in the sense of formal expansions in h (holding uniformly Vi 6 R) recalling that, in the 
foregoing section, we solved 

(0) 

H {e- Sh ' h ) = h£ n {e- Sh ' h ) (V.30) 

(0) 

for Sfr and £n in the analogous way. 



Expanding out Eq. (IV.29j) and using Eq. flV.30j) to reexpress those terms in which H acts 



purely on the e Sh ^ h factor one readily derives the associated Schrodinger equation for 0^, 

1 dS h din h d 2 %, / « \ ( V 



m dx dx 2m dx 2 



A£ n <Pk (V.31) 
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where 



(*) (*) (0) 

A£fi := £ h — £ h 



(*) 



(o) 



(*) 



(0) 



f (o) - £(<>) ) + M ^(i) - J + gf ( E m - ^(2) 



ft 2 Aw 



(I)) 



ft 3 /w 



(«) 



+ gf ( f (3) - £(3) ) + • • • 



(V.32) 



(*) (*) 
A£ (0) + hA£ (1) 



h 2 . w ft 3 . w 
¥ A^ (2 ) + ¥ A£ (3) + . . . 



I (*) (*) (*) 

designates the 'gap' between the sought-for energy coefficients < £( ), £(i), £(2), 

{(o) (0) (0) ^ 
£(0), £(1), £(2), + • • • f • 

(*) 



and those 



In the spirit of the preceeding section we now expand 0^ in the analogous fashion, setting 

(V.33) 



m (*) w fr 2 (*) h 3 (*i 

M x ) = 0(o) (x) + H(i)(x) + — (2 )(^) + gj-^(3)(a:) + • • • 



(*) 



and substitute this together with the previously defined expansions for A£ h and into 

Eq. (IV. 31 p . requiring the latter to hold, uniformly V x G M, order by order in H. 

(*) 

The transport equation for the zeroth order term 0( O ) is easily seen to be 



(*) 

1 dS {0) d0(p) 
m da; cfe 

wherein, from the preceeding section, one has 



(*) \ (*) 
A£ (0 ) 



'(0) 



mw x 4 / 1 + 



2(yfX 2 



2 • 



(V.34) 



(V.35) 



Integrating Eq. (IV. 341) one finds the general solution 



(*) 



0(0) (x) = cjo) 



1 + a/1 + 



2gx 2 



(V.36) 



where cm) is a constant of integration. This expression is everywhere smooth and nontrivial 



(*) 



if and only if C( ) 7^ and is a positive integer. Thus we henceforth refine the notation 
by writing 



AS, 



(n) 



( o) (x) := C( ) 



1 + Jl + 



(V.37) 



(n) 



Af (0) = no>o, C(o) = constant ^ 0, 
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(n) 

with n = 1, 2, 3, . . .. Note that 0(o) is bounded Vi G 1 but reduces, in the limit g \, 0, to 

(n) 

C(Q)X n which is simply a constant multiple of the leading term in the (unbounded) Hermite 
polynomial H n (\/^^ x )- Recall as well that, in this same limit, e~ Sh{ - x ^ h reduces to the 
gaussian factor e _5 l r ) x . 

It is worth nothing at this point that, had we chosen to study the (higher order) oscillators 
having potential energies given by 

TiLOnX + qx , 

(V.38) 



v *( x ) = 2 muj o x + 9 X 



k = 2,3,4, 



we would have obtained the following, more general result 



(n) 



K0(O)(z) = «C(0) 



X 



1 + 



2gx 2 ^- 1 ) 



!/(«-!) 



(V.39) 



with k = 2, 3, 4, 5, . . . and n = 1,2,3,.... Notice that, for each allowed value of k and n, 

(n) 

K 0(o) (x) is bounded V x G R but again reduces to a constant multiple of x n in the (harmonic) 
limit g \ 0. 

Returning to the quartic case and its associated Schrodinger equation ( 1V.31I) one readily 

(n) 

derives the transport equation for the first order 'correction' 0m 



1 dS, 



(n) 



(0) «</>(!) 



(n) 



m dx dx 



(n) 



(n) (n) 

A£ + 1 _ 1 dS W #(o) 



(V.40) 



wherein, of course, 0m) is a smooth, nontrivial solution to the zeroth order equation 



(n) 



1 dS {0) d<f> (0) w _ 

- nu) Q <P(o) — u, 



dx dx 



(V.41) 



(») 



of the type just discussed, S'(i) is given by Eq. (1III.7j) and A£(i) is an as yet undetermined 
constant. Comparing Eqs. (IV.40P and flV.4ip one sees that a natural technique for solving 
the first of these is the method of 'variation of parameters' wherein one seeks a solution of 
the form 



In) 



( i) (ar) = u ( i)(:r) 0(o) (x) 



(n) 



(V.42) 
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Some straightforward calculations, which made use of the preceeding, explicit expressions 

(n) 

for S'(o), S(i) and 0( O ), lead easily to the following equation for u^y. 

fr(Qfrn3 I V 9 / ^ K J J I 

(V.43) 



+ 



>Q(x)) 



w _ 2mg(nuj ) (2Q(x) + froaft) 
Q*(x)(Q(x) + mu ) 



where 



Q(x) := muQ\\l + 



mw5 



1 dS(p)(x) 
x dx 



(V.44) 



(n) 

The integral yielding can be done explicitly but generates a logarithmic singularity 

(n) 

at x = unless A£(i) is chosen to be 

v ; 2 m 2 u)Q 
in which case the expression for reduces to: 



WK J Amul{xQ{x)) 2 1 

+ gmco 2 x 2 (9n + 5n 2 ) + #V(18n + 10n 2 ) ( V - 46 ) 



— (mcjQ + 6(/a;oa; 2 )(5(a;)(n + n 2 )} . 
In deriving this we have adjusted the (otherwise arbitrary) constant of integration so as to 

(n) 

cancel a constant term in the Laurent expansion for which now takes the form 

w , , -n(n-l) 1 (25n + 9n 2 )g 2 x 2 



4muo x 2 8m 3 WQ 
{— 45n — 13n 2 )g 3 x 4 



(V.47) 





Though ti(i)(x) has, for n > 1, a singularity oc n( - n ^ at the origin, the resulting formula for 
0(i) = W(i)0( O ) is easily seen to be smooth, and in fact bounded, V x G K. In the harmonic 
limit however 

Mix) ^ -^-^^ x- 2 (V.48) 

so that 

(n) 



55 



which (for n > 1) is simply a constant multiple of the two leading terms in H n (y^par). 

Combining Eqs. f!V.37j) and flV.45j) one obtains the first order approximation for the 
energy gap 



(n) (») (0) 

AE n :— E n — E h 

(n) (n) 

~ h ( A£ (0) + hA£ {1) + ... 



(V.50) 



HuJn 



n H — 



3 gh 2 

I <n 

2 9 3 



n 2 + n) + . . . 



(n) 



The transport equation for the second order 'correction' 0( 2 )(x) is readily computed to 



be 



1 dS, 



(n) 



(0) a<P(2) 



(n) 



m cfe dx 



nwo0(2) 



(n) (n) (n) 

1 dS(2) d<p(o) 2 (iS'(i) d(j)(i) 1 d 0(i) 



+ 



(V.51) 



m dx dx m dx dx m dx 



(n) (n) (n) 

+ A£( 2 )0(o) + 20(i) 



3 gn(n + 1) 
2 m 2 WQ 



wherein 



(n) 



(») (n) 

0(o), 0(i), 5(0), 5*(i) f are the functions defined above, £(2) is given by Eq. (1V.13j) 



and A£p) is an an yet undetermined constant. 

Applying the method of variation of paramters we set 



(n) 



(n) 

0(2) (X) = W(2)(X)0(Q)(Z) 



(n) 



(V.52) 



'"(2) 



(n) 



and derive, easily, the associated explicit formula for —4-*-. The integral yielding U( 2 ) contains 



(n) 



a logarithmic singularity at x = unless A£( 2 ) is chosen to be 



(59n + 51n 2 + Un 6 ) 



(V.53) 
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in which case the expression for um) reduces to: 

+ 3gm 3 u%x 2 (35 + 28to + 32to 2 + 5to 3 ) 

+ 2g 2 m 2 cu*x\-257 - A2n + 104to 2 + 75n 3 ) 

+ 28# 3 mu; 2 a; 6 (-59 - 24to + 8n 2 + 15n 3 ) 



+ 24C/VX-59 - 24to + 8to z + 15n d yj <y . j , 

+ Q{x) [3m 4 ^(16 + 2 In + 2to 2 + to 3 ) 
+ Qgm 3 ujQX 2 (12 + 37to + 24to 2 + In 3 ) 
+ # 2 m 2 c4r 4 (-1211 - 351to + 380to 2 + 282to 3 ) 
+ 4# 3 mw 2 x 6 (-1355 - 837to + 8to 2 + 156to 3 ) 
+ 4# 4 x 8 (-1355 - 891to - IOOto 2 + 102to 3 )] } 

wherein, as above Q(x) is given by Eq. (IV.44j) . We have again adjusted the choice of a 



(n) 

constant of integration so as to cancel a constant term in the Laurent expansion for M( 2 ) 
which now takes the form 

(n) (to — 3)(n — 2) (to — l)n 1 



TO (2 )(XJ 



16toi 2 o;q x 4 



g to 2 (to - 1) 

m 3 u; 4 x 2 
3 

+ ~ t « (-1098 - 719to - 140to 2 + 13to 3 )o; 2 ( V - 55 ) 
16tto 5 Wq 
4 

(7 TO 

+ - fi 10 (12083 + 7217to + 1460to 2 - 4to 3 )x 4 



64m 6 



+ 0[x 6 ]. 

Since the singular term oc ( n ~ 3 )( n ~^)( w ~ 1 ) w j s present only for to > 4 and that oc - ^" 2 ~ 1 ^ only 

(n) (n) 

when to > 2 it's clear that the product (f)(2) = TO(2)</>(o) is m fact smooth V x G K. It is easily 
seen to be bounded as well though, in the harmonic limit it reduces to 



\ ( n ) 

- 2) (x) — > M"-3)(n-2)(n-l) ^ 
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so that 



( ? J? h2(n J \, ^ C(o) f n n(n-l) / ft \ n _ 2 

n(n- l)(n-2)(n-3) / ft \ 2 
32 \mw / 



(V.57) 



n— 4 



which (for n > 3) is a constant multiple of the three leading terms in H n (v/^Tp^) • ^ is not 
difficult to prove that, in this limiting case (i.e., setting g = everywhere), our method will 
simply regenerate the full Hermite polynomials and then terminate thereby reconstructing 
the exact excited states (after multiplication by the gaussian factor) for the pure harmonic 
oscillator. 

Returning to the quartic oscillator results and combining Eqs. (1V.37I) . (1V.45I) and (TV. 531) 
we find that our approximation yields the following second order formula for the n-th excited 
state's (total) energy 



(n) 



1\ 3 / gh \ f 2 1 
n + - + - -Ar^ [n 2 + n + - 
2/ 2 \m 2 ujl) V 2 



I ( gh ^ ' 



(34n 3 + 51n 2 + 59n + 21) 



(V.58) 



This agrees with the standard result computed using conventional perturbation theory car- 
ried out to the corresponding order [yj]. We conjecture that such agreement will persist to 
arbitrary order but a full proof of this, if true, would presumably require the development of 
suitable recurrence relations for the (excited state) energy coefficients. We shall not pursue 
that issue herein. 

It may seem paradoxical that, to the orders that we have computed (and conjecturally to 
all orders), both ground and excited state energy formulas agree with those of conventional 
perturbation theory whereas the corresponding wave functions of the two approaches differ 
dramatically. Our expectation though is that one could use our results to regenerate the 
conventional approximate wave functions of perturbation theory by expanding our formulas 
for Sfr(x), etc., in powers of g, Taylor expanding the coefficients of the gaussian factor 
e ~2\~h~) x a nd truncating the results at the desired order in g and finally reexpressing these 
truncated coefficients as finite series of Hermite polynomials |iif n (a/^tpx) |. 

(n) 

We conclude this section by pointing out that Eq. (1V.39j) for re 0( O )(x) simply expresses 
this quantity as (a constant multiple of) the n-th power of the corresponding Sternberg 
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coordinate for the associated oscillator. More precisely, the transformation 

2 i/(*-i) a . 



(1 + /i + i/( w -i) 



with smooth inverse 



maps l9i diffeomorphically to the interval 



and transforms the Hamilton- Jacobi flow vector field 



(V.59) 



X fl a v 2( k -i))i/(k-i) ( V - 60 ) 



3 y (V.61) 



X H := -^-=w xi H 2—~T (V.62) 

m dx dx y mwg dx 

to the Sternberg form 

X H = oo y^- (V.63) 
"1/ 

that was discussed, in a more general setting, at the end of Sect. HVl 

Though we did not exploit this transformation to carry out the foregoing calculations it 
is quite conceivable that they would have been simplified thereby and, furthermore, that the 
use of Sternberg coordinates could facilitate the resolution of some of the conjectural issues 
mentioned above. 



VI. CONCLUDING REMARKS 



For the Lagrangians normally considered in classical mechanics it would not be feasible to 
define their corresponding action functionals over (semi-) infinite domains, as we have done, 
since the integrals involved, when evaluated on solutions to the Euler Lagrange equations, 
would almost never converge. It is only because of the special nature of our problem, with 
its inverted potential energy function and associated boundary conditions, that we could 
define a convergent action integral for the class of curves of interest and use this functional 
to determine corresponding minimizers. Since this problem is of a highly non-standard 
type within the calculus of variations we have felt obligated to give, in Sect. a rather 

complete and self-contained analysis of the existence, uniqueness and smoothness properties 
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of the associated minimizing curves together with a proof of smoothness and asymptotic 
behavior for the action function, SVo) (x), computed from these minimizers. 

A remarkable feature of our construction, given the hypotheses of convexity and coercivity 
imposed upon the potential energy V(x), was that it led to a globally smooth solution to the 
corresponding ipze Hamilton- Jacobi equation. Normally the solutions to a Hamilton- Jacobi 
equation in mechanics fail to exist globally, even for rather elementary problems, because 
of the occurrence of caustics in the associated families of solution curves. For our problem 
however caustics were non-existent for the (semi-)flow generated by the gradient of S , ( )(x). 
The basic reason for this was the inverted potential character of the forces considered which 
led to the development of diverging (in the future time direction) solution curves having, 
in effect, uniformly positive Lyapunov exponents that served to prevent the occurrence of 
caustics altogether. 

By contrast, the more conventional approach (in the physics literature) to semi-classical 
methods leads instead to the standard (non-inverted-potential-non-zero-energy) Hamilton- 
Jacobi equation for which, especially in higher dimensions, caustics are virtually unavoidable 
and for which, even in their absence, a non-trivial matching of solutions across the boundary 
separating classically allowed and classically forbidden regions must be performed. While 
Maslov and others have developed elegant methods for dealing with these complications 
[lsl ] their techniques and results are more appropriate in the short wavelength limit wherein 
wave packets of highly excited states are the central objects of interest. On the other hand 
our approach is aimed at the ground and lower excited states though, in principle, it is not 
limited thereto. 

As we have already mentioned though, our approach is a natural variation of one that 
has been extensively developed in the microlocal analysis literature but it also differs from 
this innovative work in fundamental ways that are crucial for our ultimate, intended appli- 
cations. In the microlocal approach [2|, one begins by analyzing the (classical, inverted 



potential) dynamics 
rem of mechanics 



ocally, near an equilibrium, by appealing to the stable manifold theo- 
16| . One then shows, by a separate argument, that, for an equilibrium p 
(lying in some neighborhood U C R n ) the corresponding stable (W s (p) C T*U) and unsta- 
ble (W u (p) C T*U) submanifolds of the associated phase space T*U 8X6 ill fcLCt LcLgrcLllgicLll 
submanifolds that can be characterized as the graphs of the (positive and negative) gradients 
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of a smooth function : U — > R: 

W U { P ) = {(x, p) | x e U, P = V0(x)} 

= {(x, p) I x £ Z7, p = -V0(x)} . 

This function is shown to satisfy a certain 'eikonal' equation (equivalent to our ipze 
Hamilton- Jacobi equation restricted to U) and 0(x) itself is, of course, nothing but the 
(locally defined) analogue of our action function 5 , ( )(x). 

The potential energies, V(x), dealt with in the microlocal literature often entail multiple 
local minima, or 'wells', for which our global convexity and coercivity hypotheses are not 
appropriate. Much of the detailed analysis therein involves a careful matching of locally 
defined approximate solutions (constructed on suitable neighborhoods of each well) to yield 
global asymptotic approximations to the eigenvalues and eigenf unctions for such problems. 
Since, however, we are focused primarily on potential energies having single wells (corre- 
sponding to unique classical 'ground states'), many of the technical features of this elegant 
analysis are not directly relevant to the issues dealt with herein. 

For this case of a single well, however, we have essentially unified and globalized several 
of the aforementioned, local arguments, replacing them with an integrated study of the 
properties of the (inverted potential) action functional T ip [y} {cj., Eq. flHI.lj) ). When one 
turns from finite dimensional problems to field theoretic ones |l| this change of analytical 
strategy begins to play a crucial role. For the typical (relativistic, bosonic) field theories of 
interest to us in this context, the Euler Lagrange equations for the corresponding, inverted- 
potential action functionals that arise are the Euclidean signature, elliptic analogues of the 
Lorentzian signature, hyperbolic field equations that one is endeavoring to quantize. While 
generalizations of the aforementioned stable manifold theorem do exist for certain types of 
infinite dimensional dynamical systems, the elliptic field equations of interest to us do not 
correspond to well-defined dynamical systems at all. In particular their associated Canchy 
initial value problems are never well-posed. 

On the other hand certain elliptic boundary value problems for such Euclidean signature 
action functionals are mathematically meaningful and, fortuitousl y, a number of the most 
important of these have already been the object of rigorous study Q]- In a companion 
paper to the present one we shall review some of the principal results of this analysis and 
extend them with the aim of generalizing the arguments given herein to a natural infinite 
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dimensional setting 
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Appendix A: Appendix 

Sternberg's linearization theorem jsj implies that, at least on a neighborhood U of the 
origin in M n , there exists a diffeomorphism 

x^ A1 (x) = (y 1 (x),...,y»(x)) (A.1) 



with 



and 



y\*)=x i + 0{\x\ 2 ) (A.2) 
dy l 



v S 

such that on n(U) the vector field (0) takes the form 



5) + 0{\x\) (A.3) 

V5„ 



y j.ayx) a _^ y , a 

z — ' m dx % dx % t* ^— ' dy % 

i=i i=i y 

By exploiting special features already established for (0) however, we shall be able to 
extend Sternberg coordinates to a global chart defined on a star-shaped domain K C K" 
with iT x (K) « M n = {(x 1 , . . . , x n )}, such that 



I VS(0) = I>^ (A.5) 

1=1 " 

everywhere on K. For this purpose we shall first need some further insight into the (semi-) 
flow generated by -jrVS^o). 
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1. Global Properties of the Gradient (Semi-) Flow of Sr \ 

We know from Sect. (1IIII) that SVq) : M. n — > R is a globally smooth solution to the ipze 
Hamilton- Jacobi equation 

1 ^fdS {0) (x) 



i=l 

with the asymptotic behavior 



5 (0) (x) = im^a; i (x0 2 + O(|x 



i=l 



(A.7) 



d^o^x) = mujjX 3 +0(\x 



2\ 



This function was determined by minimizing the ip action functional (lIH.lj) on the affine 



space of curves £> x (c.f., (IIII.2I) ) having (arbitrary) fixed endpoint x e M n , evaluating the 
action on the minimizing curve, 7 X , and setting 

S (0) (x) :=X ip [ 7 x]. (A.8) 

In view of the coercivity condition ( 1III.20I) assumed for the potential energy however, it 
follows that, for any curve 7 (t) = (x x (t), . . . , x n (t)) lying in D x , one has 

T tp [l] > VM (A.9) 



where 



with 



X;[ 7 ]:= J |^m|j{(^(t)) 2 + ^(^(t)) 2 }|^ (A.10) 



: \Iuj1 - A, 2 > V i E [1, . . . , n]. (A.ll) 
But the minimizing curves for X* p within D x are simply given by 

7 *(t) = (xV 1 *, . . . , x n e Vnt ) (A. 12) 

for which one finds that 

S(* )(x) :=V P K] 

1 " 
= -m^Uiix 1 ) 2 . 



(A. 13) 
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It follows that, globally on 



<S(o)(x) > 5* 0) (x) 



~mf>(^) 2 (A.14) 



2 



and thus, in particular, that any level set of with level value less than s > lies inside 
the (hyper-) ellipsoid defined by 

S?o)(*) = *. (A.15) 

From Eq. ( IH.lal) and the properties assumed for V(x) it's clear that the gradient of S(q) 
only vanishes at the origin which coincides with the (exceptional) level set 5^(0). For any 
level value s > the corresponding level set, Sz}(s), is therefore necessarily a smoothly 
embedded, n — 1 dimensional submanifold of IR n . In view of its boundedness, Szl(s 



is 



compact and, from the asymptotic behavior of near the origin (c.f. Eq. (1III.68j) ), it's 
clear that the level sets for sufficiently small s > are topological n — 1 spheres surrounding 
the origin. Since V5*(o) is nowhere vanishing away from the origin it moreover follows, from 
basic Morse theory, that every level set of S'(o) (corresponding to a level value > 0) is also a 
smoothly embedded, topological n — 1 sphere. 

In view of the fact that, along the integral curves 7 X of ^VS'(o), £(0) obeys 

4s ( o) (7x(t)) = -J-VS ( o) ■ V5 (0 ) ( 7x (i)) 

dt m (A. 16) 

= 2F(Tx(*)) 

we see, since V(x) > except at the origin, that S'(o) (7x(£)) is monotonically decreasing as 
t \ — oo with 

lim S (0) ( 7x (t)) = % (0, . . . , 0) = VxGl". (A.17) 

t\— oo 

Let B s C M. n designate, for any s > 0, the open n-ball bounded by the (topological) 
(n— l)-sphere S7 Q As). From the asymptotic behavior of 5*(o)(x) (c.f. Eq. flA.7j) ) there clearly 
exists an s* > 0, sufficiently small, such that B s * lies entirely within the neighborhood 
U C R n on which Sternberg's theorem ensures the existence of the diffemorphism \x described 
above. From the monotonicity property (1A.16l4A~.17p every integral curve of — V5(o), followed 
in the negative time direction, eventually enters and remains within B s *, asymptotically 
approaching the origin as t \ — oo. 

In fact, using the monotonicity formula, it is easy to place an upper bound on the time 
required for any point x G B s to flow to (and remains inside) the 'Stenberg domain' B s *. For 
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s < s* there is nothing to do so assume s > s* and let A s * yS C M™ designate the compact, 
annular solid region defined by 

A s *,s := K\B S * (A.18) 

and thus bounded inside and out by S^(s*) and S^(s) respectively 

Since the potential energy function V(x) is continuous and strictly positive on the com- 
pact annulus A s *, s it achieves (strictly positive) maximal and minimal values thereon. Let 
V™ 1 ^ > designate the corresponding minimum. From (1A. 16|) it follows that, for any 



t < := 2 y m f n < 0, every integral curve 7 having 7(0) G B s will satisfy j(t) G B s *. Thus, 

s* ,s 

along the (semi-) flow generated by — VS(o), B s gets mapped to an image lying within B s * 
for all t < T™™ < 0. 

By standard results on the (semi-) flows of smooth vector fields whose integral curves are 
all complete in the negative t-direction, it follows that the aforementioned mapping is, in 
fact, a C°° diffeomorphism [191 ] . Thus 



7* : B s — > 7*(B a ) C B„. (A.19) 



defined, for t < T^ a s x , by 



7*(x)= 7 x(*), (A.20) 

where 7 X ( - ) is the integral curve of — ViS'(o) satisfying 7 X (0) = x G B s , is a diffeomorphism 
taking B s to its range lying within the Sternberg domain B s *. 



2. Extending the Sternberg Diffeomorphism 

Let 7x(^i), 7x(*2) designate any two points along an arbitrary integral curve 7 X of — V^o) 
that both lie within the Sternberg domain B s * and assume t\ < t?, < 0. Note that t2 = 
would be allowed here if x = 7 X (0) already lies within B s * . Utilizing the Sternberg coordinate 
expression for ^■VS'(o), given by Eq. ( 1A.4j) and valid throughout B s *, one easily shows that 

// (7x(ti)) e"^ 1 = // ( 7x (t 2 )) e"^* 2 (A.21) 

V i G [1, . . . , n] and for any such choice of 7 X , t\ and t 2 - 

Let us tentatively extend the domain of definition of the Sternsberg map fi (assumed for 
simplicity to coincide initially with the n-ball B s *) to B s by setting, for any x G B s and for 
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some t < T™ x < 0, 



y l (x) = //(x):=//( 7x (t))e-^ (A.22) 



V i G [l,...,n]. In view of Eq. ( 1A.21I) it is easily verified that the map so defined is 
independent of the actual value of t chosen. Furthermore since, for any such t, 7*(-) = 7-(i) 
is a diffeomorphism defined on B s , we see that the formula 

y\.)=^{.):=^{ 1 \.))e~^ t (A.23) 

V i G [1, . . . ,n], for the extension of fi to B s , is simply a composition of diffeomorphisms. 
By virtue of flA.211) it is easily checked that this extended fi, if restricted to B s *, coincides 



with the Sternberg map originally given. 

One can repeat this construction for a sequence Sj — > oo whose corresponding domains 
B Si exhaust M. n as i — > oo and verify as above that, with each successive enlargement of the 
domain for /i, its restriction to the previous domain coincides with the previous definition. 
In this way one arrives at a resulting diffeomorphism fi whose maximal domain of definition 
exhausts M n . 

To see that this extended /i map has the desired property of transforming ^VS'(o) to the 
Sternberg form (1A.5I) consider an arbitrary segment 

7x : (— e, e) — > (A.24) 

of an arbitrary integral curve 7 X of — V^o). For any e > one can choose a T < sufficiently 
large and negative that 

7T*(A)(T):(-e,e)— >R» (A.25) 

will be entirely within B s *. By virtue of the autonomous character of the vector field ^VS'(o) 
it follows that 

7 7x (A)(T) =7 7x{0 )(T + A) 

= l x (T + \)eB s * 
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V A G (— e, e). Computing the Sternberg coordinate form for this curve segment we get 

y%\) :=^( 7x (A)) = /i' ( 77x{A) (r))e-^ T 
= // ( 7x (T + A)K- T 

= ft (7x(T + A)) e ^(T+x) e ^x ( A .27) 

= // ( 7x (T))e-^ T e^ A 

= //(x)e WlA Vie[l,...,n] 

Thus 

y*(A) := ^ ( 7x (A)) 
= //(x)e^ A 

which gives immediately that 



(A.28) 



along the arbitrary curve segment. Since /i is a globally defined diffeomorphism on lR n it 
follows that it transforms — VS(o) to the Sternberg form (1A.5[) . From the standard definition 
(c.f., p. 23 of both M™ and its image under [i are star-shaped domains for the (semi-) 
flow generated by ^VS'(o). 

We conclude this appendix with a brief discussion of how the Jacobian determinant of the 
transformation to Sternberg coordinates varies along the integral curves of — VSVo). Writing 
y l = [^(x 1 , . . . , x 11 ) and {x a } = {x 1 , . . . , x n } we have, 

1 detfM) (A.30) 



where 

is the usual transformation relating the Euclidean metric e in Sternberg coordinates, 

e = gij (y) dy l <g> dy j (A.32) 

to its Cartesian form 

e = 8 ab dx a <g> dx b . (A.33) 
A straightforward calculation results in 



^ln(v/d^) = |lndet (g^) 



(A.34) 

bji. 



m 

i=i 
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